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Abstract. We study a low temperature anisotropic anti-ferromagnetic 2D Ising model 
through the guise of a certain dimer model. This model has a bijection with a one- 
dimensional particle system equipped with creation and annihilation. In the thermody- 
namic limit, we determine the explicit phase diagrams as functions of temperature and 
anisotropy. Two values of the anisotropy are of particular interest - the 'critical' value and 
the 'independent' value. At independence, the particle system has the same distribution as 
the two colored noisy voter model. Its limiting measure under a natural scaling window is 
the Continuum Noisy Voter Model. At criticality, the distribution of particles on a given 
horizontal line, is a Pfafnan point process whose kernel in the scaling window can be written 
explicitly in terms of Bessel functions. We also show that 'macroscopic' creations form a 
Poisson point process. 



1. Introduction 

1.1. History. The dimer model is a two-dimensional exactly solvable model, which can be 
described in the following manner. A dimer configuration for a planar graph is a subset 
of edges in which each vertex is covered exactly once by an edge. The dimer model is a 
probability measure on the set of all dimer configurations of the underlying planar graph. 
Weights can be applied to the edges so that specific edges have higher or lower probabilities 
of being observed. 

The dimer model was initially introduced in [Kas61j. [Fis66j noticed the correspondence 
between the dimer model and the two-dimensional Ising Model. Since its introduction, 
the dimer model has been studied extensively and a vast variety of techniques have been 
developed using many areas of mathematics. Much of this work was completed by R. Kenyon 
and A. Okounkov during the late 90's and 00's (for an excellent survey of their body of results, 
see |Ken09] ) . The dimer model is a rather unique statistical mechanical model - not only 
can a partition function be easily found but by [Ken97| . local statistics can be computed. 

Dimer models on non-bipartite lattices are fundamentally different to dimer models on 
bipartite lattices. For instance, we can define a height function if and only if the underlying 
lattice is bipartite (see [KenOO]), which leads to connections with random surfaces ([Shc05]). 
Recently, there has been renewed interest in the study of dimer coverings on non-bipartite 
lattices. In particular, [BdT08j have shown that some of the machinery for the dimer model 
on bipartite lattices transfers to non-bipartite lattices. 

1.2. The Dimer model on the Fisher lattice. This paper focuses on a dimer model on a 
particular non-bipartite lattice, called the Fisher Lattice, introduced in [Fis66] . The Fisher 
lattice consists of a honeycomb lattice with a triangle decoration at each vertex as depicted 
in Figure [TJ All vertical edges are referred to as a edges. All edges outside the triangle 
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Figure 1: The Fisher Lattice 



decoration which are not vertical are referred to as b edges. We give weight a to the a edges, 
weight b to the b edges and weight 1 to the remaining edges. 

Let G denote the Fisher lattice on the plane. Define the toroidal graph G n to be the 
quotient Gjnl? where the action nZ 2 are translations by (n,0) and (0, n) (horizontal and 
vertical translations). Let M{G n ) denote the set of all dimer coverings of G n . We can define 
a probability measure fi n for M G M(G n ) so that the probability of M is proportional to 
the product of all the edge weights; that is, 

a N a uN b 

(1.1) ,„(M) = ^ 

where N a = N a {M) and = Nb(M) are the number of dimers covering a and b edges for 
the dimer covering M and, Z(G n ) is the partition function: Z(G n ) = X^A/eX(G„) aNa b Nb - 



1.3. Anti-ferromagnetic Ising model and Setup. Here, we descibe the one-to-two corre- 
spondence between the dimer model on the Fisher lattice and an anisotropic anti- ferromagnetic 
2D Ising model (explained below), as described in [Fis66j. 
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Let A n represent the faces of the dodecahedrons in G n . Let £1 := {— 1,+1} A ™ and let 
a x G {—1, 1} denote the spin at x G A n . For a € f2, define the Hamiltonian by 

(1.2) ^) = EE J ^^i ±1 

where u ~ w means that v and w are nearest neighbors and J vw is the interaction term 
between the faces i> and w. For d£!l, define the probability measure, 

e -H(a) 

(1.3) ,„M = ^ 

where Z(A n ) is the partition function: Z(A n ) = Yla^n e ~ H ■ The a bove probability mea- 
sure is called the Ising model. 

The correspondence between the Ising model and dimer model can be seen in Figure [2] 
and is as follows. Faces sharing an edge have the same spin if there is a dimer covering the 
shared edge. Otherwise, the faces have opposite spins. 

log a w ~ v and wv is an a edge 
log b w ~ v and wv is a b edge, 



(1.4) J v 



where wv is the shared edge (on the Fisher lattice) between the faces w and v. It is clear 
that 

Lemma 1. 

(1.5) 2Z(G n ) = Z(K n ) 

We choose a = x and b = ux where x < 1 and u £ (0, 1). As u < 1/x then J VjW < for 
v ~ w, which means that adjacent sites have a high probability of having opposite spins. This 
is called an anti-ferromagnetic interaction. Furthermore, J VjW depends on the direction of the 
edge between the sites v and w. This is called an anisotropic interaction. The parameter x 
controls the temperature (because the temperature is equal to l/(log(l/x))). The anisotropy 
is controlled by the parameter u. We shall sometimes refer to u as the anisotropy. 

Finally, we heuristically describe the concept of thermodynamic limit. As previously men- 
tioned, dimer coverings of the torus G n give a probability measure fi n . The thermodynamic 
limit is the measure, fj, = fj,(u,x) = linin-^oo /i n , i.e. the limiting measure obtained when the 
system size is sent to infinity. More information on the thermodynamic limit can be found 
in Section [2) 

1.4. The particle model. By considering the horizontal axis as space and the vertical axis 
as time, the dimer model, fi, is in bijection with a one-dimensional particle system called the 
particle model. This is a probability measure on the set of all possible bi-infinite lines, finite 
loops and pairs of infinite lines with the same starting points on the diagonal grid. 

For each dimer configuration there is a realization of the particle model, and vice versa, 
chosen with the same probability. An example of the correspondence can be seen in Figure 
El Each realization of the particle model is obtained from a dimer covering of the Fisher 
lattice by collapsing the a edges and the triangle decoration. This leaves a diagonal grid, 
where each dimer covering a b edge in the dimer model can be thought of as a trajectory of 
a particle in the particle model. By the structure of the underlying dimer model, particles 
can be created or annihilated in pairs and there is no branching or coalescence of particles. 
The above construction shows that the particle model is in bijection with dimer model of 
the Fisher lattice. 
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Figure 2: The correspondence between the Ising model and dimer coverings on the Fisher 
lattice. The blue edges represent dimers covering b edges. The red edges are dimers covering 
the decoration and a edges. The '+' and '-' are spins in the corresponding Ising model. 



We say that a singularity occurs when there are no dimers covering any edges at a triangle 
decoration. A singularity is called a creation ( annihilation) if the a edge is below (above) the 
triangle decoration. We use the terminology creation (annihilation) because that singularity 
is a creation (annihilation) of two particles in the particle model. We adopt the convention 
that a particle is located at a point if there is no dimer covering the corresponding a edge. 

1.5. Scaling Limits and Scaling Windows. Here, we heuristically describe the concepts 
of scaling limit and scaling window. More specific information regarding scaling windows for 
the dimer model on the Fisher lattice, can be found in Section [21 

The scaling limit of grid-based models is a subtle concept and we make no attempt to 
define it rigorously. Informally, we can think of the scaling limit as the limiting measure 
of some measure when sending the lattice spacing to zero. For a two-dimensional lattice 
(e.g. Z 2 ), the vertical and horizontal lattice spacings need not be sent to zero at the same 
rate. For example, Donsker's scaling of a simple symmetric one dimensional random walk 
has different vertical and horizontal scalings (by considering the trajectory of the walk on a 
diagonal grid). For many statistical models, scaling limits are only non-trivial provided that 
the model is at criticality. In this context, criticality refers to the correlation length, defined 
in (j!.16j) . being infinite. See [FFS92] for discussions regarding criticality and [SchOO] for the 
underlying philosophy of scaling limits of statistical mechanical models. 
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Figure 3: A dimer configuration with its corresponding particle model 



The concept of the scaling window is equally as subtle as the concept of the scaling limit. 
In this case, we do provide a watertight mathematical definition for the particle model (see 
Section [2|) as it is one of the main aspects of this paper. Roughly speaking, a scaling window 
is the limiting measure of some measure when sending the lattice spacing to zero while 
simultaneously sending a certain parameter to zero. Similar to a scaling limit, the choice of 
the re-scaling is important in guaranteeing non-trivial behavior. This choice usually involves 
having constant density of some observable in the scaling window. For example, we could 
choose the scaling window so that the density of particles is constant along any arbitrary 
horizontal line. 

Contrary to the scaling limit, the scaling window preserves microscopic behavior of the 
discrete model, which means that the scaling window is non-trivial regardless of whether the 
model is critical. In particular, the scaling window of a critical model is still critical. Indeed, 
the scaling window of a critical model can only be achieved if the model has two (or more) 
parameters. This is a consequence of criticality being a relationship between the available 
parameters. 

In this paper, all scaling windows are taken with respect to the parameter x. We will 
only consider scaling windows after the thermodynamic limit has been taken. We will only 
consider scaling windows of the particle model and the distributions of the particles on an 
arbitrary horizontal line. We use the following terminology for scaling windows: For a, (3 > 0, 
the scaling window (x a ,x^) of a measure is the limiting measure obtained when sending a 
box of lattice points of size (l/x a , 1/x 13 ) to a box in the continuum of size (1,1), when x — > 0. 

We can now state and explain our main results. First, we consider the behavior in the 
thermodynamic limit and in particular, highlight the behavior at low temperature (i.e. which 
means x is small) for different values of the anisotropy, u. We can then focus on taking scaling 
windows for different values of the anisotropy parameter. 
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1 .6. Statement of the results. 
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Figure 4: The Phase Diagram in terms of u and x. The solid line represents criticality, 
where the partition function is non-analytic. The dashed line represents independence. 



= [i(x,u) is the thermodynamic limit. Define u c := (— 1 + \/l + x 2 )/x 2 , 
x 2 )/x 2 , fio = lim^o fi(u,x) and P At (-) denote the probability of the event • 



Recall that /j, 
Ui := (1 - VT^ 
under /x. 

Under /i, Theorem [5] gives an explicit formula for the expected number of a, b edges and 
creations per fundamental domain for x 6 (0, 1) and u G (0, 1). Essentially, each expectation 
can be written as a linear combination of elliptic integrals, with the exact expression depen- 
dent on whether u < u c , u = u c or u > u c . By considering a low temperature expansion, 
then 



Theorem 1. Under [Iq, 

(1.6) P ^(particle) = 
and 

(1.7) P ^(creation) 



arccos 



if u < | 

("£) ^ U> \ 



0. 



In this limit, there is a phase transition at u = 1/2. As there are no creations (or 
annihilations), then particles form non-intersecting lattice paths for u > 1/2. The above 
behavior is only exhibited under jiq,. 

We can now consider the phase behavior for u < u c chosen independently of x which leads 
to choosing a natural re-scaling for the scaling window of the particle model: 

Theorem 2. For u < u c chosen independently ofx, the scaling window (x,x) of the particle 
model is a Poisson Point process with intensity u 2 ( l/y/T— 4u 2 — 1 ) . 



Under [i = fi(ui,x) with x > 0, the locations of particles in the particle model along a 
horizontal line are distributed according to an i.i.d. Bernoulli with probability x/(l + x). 
Each trajectory of a particle is independent of the trajectories of other particles and along 
each horizontal line, creations are i.i.d. Bernoulli with probability cx 2 (where c is some 
constant). This leads to the noisy voter model interpretation of the particle model which is 
described below. 
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The (non-noisy) voter model on Z with two colors is a time dependent probability measure 
with state space {R,B} Z . The update of each site is given by randomly choosing the color 
of a neighboring site (see [Lig85 1 ) . The noisy voter model of two colors with noise p was 
introduced in }GM95| . Each site chooses at random the color of a neighboring site with 
probability 1 — p or flips its color with probability p. [FINR06J constructed the scaling 
window (x,x 2 ) of the noisy voter model with noise cx 2 (c is a constant) and called it the 
Continuum Noisy Voter Model (CNVM). This construction relies on using the Brownian Web 
(sec [FINR04J) and the dual Brownian Web. The correspondence for the particle model at 
the special case when u = U{ and the noisy voter model is as follows. The paths in the 
particle model are shown, in Section [6J to represent the boundaries between the two colors 
in the noisy voter model. It is not surprising that 

Theorem 3. At u = Ui, the particle model has the same distribution as the color boundaries 
of the two color noisy voter model. In the scaling window (x, x 2 ), the particle model converges 
in distribution to the CNVM. 



We can now consider the behavior for n 7 := (1 — yl — x 2 ^/) / {^x 2 ) = 1/2 + ^/x 2 /8 + 0(x ). 
Notice that U-i = u c and u\ = u^. We find the stationary measure for particle locations at 
time t £ 1 in the scaling window (x, x 2 ). 

Theorem 4. For u = u 7 but 7 ^ 1, in the scaling window (x,x 2 ) the distributions of the 
particles along an arbitrary horizontal line are given by a Pfaffian point process,^, with 
kernel 



(1-8) M%h,v) 



e( 7 ) 



j,x Kyi y j y _ e ( 7 ) o 

for y E M and for y\ S M and for y 2 £ R with y\ < y2 

-Ej( yi -y 2 \) -Ej(\ yi -y 2 \ 



(1.9) MZ%( yi ,y 2 



where for ei ( 7 ) = 2 + ^2(1- 7), e 2 ( 7 ) = (2 - ^2(1 - t))I 7 >-i + (-2 + V2(l - 7))I 7 <-i 
and a > 0, we have 

(mo) e M _ i r w ^±£ 

™Jea(>y) 2^4 + 87 + 4 7 2 - 12p 2 + 4jp 2 + p 4 

1 r e i(l) 2ne~ a P 
(1.11) Ej(a) = -- / PC =dp 

" ie 2 ( 7 ) V4 + 87 + 4 7 2 - 12p 2 + 4 7 p 2 + p 4 

and 

Yy/ a _ 1 r i(7) (-2-2 7 -p 2 )e- a P 



(1.12) El (a) = — / I g y dp. 

«ie 2 ( 7 ) 2 ^4 + 87 + 4 7 2 - 12p 2 + 4 7 p 2 + p 4 

In i/ie special case when 7 = —1, e(— 1) = 2/7T, 

(1.13) ^"V) = ^/(0,4a) - S L (0,4a) = - T e" 4Qsine ^ 
and 

(1.14) ^ 2 _1 (a) = 5 L (-l,4a) - 5/(0, 4a) = — C e ie - Aasine d6 

™ Jo 
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Figure 5: A plot of C'(0.2, 7 ), defined in (|1.15|), for g £ (-2,2). 

where Bi(n,z) is the modified Bessel Function of the first kind of order n and Sl(i%,z) is 
the modified Struve function of order n. 



A Pfaffian point process is the analog of a determinantal point process for an anti- 
symmetric matrix whose entries are given by some integral kernel, with the determinant 
replaced by a PfafHan. More information about determinantal processes can be found in 
[SobOOJ . 

The above theorem gives some information about the covariance between particles located 
on a horizontal line. From Theorem HI the covariance between two particles on the same 
horizontal line is 

(1.15) C(a, 7) := Ej (a) 2 - Ej(a) 2 

for 7^1, where a is the distance between two particles. For 7 = 1, C(a, 1) = for all 
a > 0. Figure [5] shows a plot of this function for a = 0.2. It can be easily shown that C(a, 7) 
is positive if 7 < 1, is negative if 7 > 1 and if 7 = 1. 
As in [CCFS89| . define the correlat ion length to be 

(1.16) £(7) == " lim r—fr^ rrr. 

a^oo log(|C(a,7)|) 

This measures the rate of decay of the covariance. In Section we prove the following 
lemma. 



Lemma 2. 



(1.17) 



£(7) 



1 

2e 2 (7) 
OO 

1 

2e 2 (7) 


1 
1 



7 < -1 
7 = -l 
-1 < 7 < 1 
7 = 1 
7 > 1 



where 62(7) is defined in Theorem^ 

The discontinuity at 7 = 1 is due to the definition of the correlation length as opposed to 
any explicit change in phase behavior. This is explained in more detail in Section [71 One 
further interesting consequence of the above lemma is the following: We have written the 
model at criticality (i.e. 7 = — 1) in terms of the parameter x and taken the scaling window 
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(x,x 2 ). This is an example of the scaling window preserving scale invariance. Here, scale 
invariance means that the correlation length diverges. 

For u = u-y, the creations (and annihilations) forms a dense set regardless of the value of 
7 in the scaling window (x,x 2 ). This is a consequence of 

(1.18) of Creations per fundamental domain] = Q(x 2 ) 

for fi = /i(u 7 ,x). In particular, for 7 = — 1 (or u = u c ), define T e to be the set of creations 
which give paths of length greater than e in the scaling window (x,x 2 ). Provided e > 0, we 
show that T e is locally finite. In Theorem [6j we show that T e is a Poisson point process with 
an intensity dependent on e. Unfortunately, our method does not give a specific formula for 
the intensity. 

We have determined the scaling window (x, x 2 ) of the particle model for 7 = 1. We have 
also found the distributions of particles along a horizontal line in the scaling window (x, x 2 ) 
for all 7 G R. We ask the following question, for each « 7 , 7 / 1, what is the scaling window 
(x,x 2 ) of the particle model? The difficulty lies in the fact that the underlying trajectories 
of particles are no longer independent, whereas the underlying CNVM construction relies 
on coalescing but independent random walks (the random walks are independent up to 
coalescence) . 

1.7. Overview of the paper. We begin by giving a more rigorous overview of the dimer 
model on the Fisher lattice, thermodynamic limit and scaling windows. We also provide 
some calculations concerning Pfaffians that are useful later in the paper. This is all con- 
tained in Section [2j In Section [3j we calculate exactly the expected number of a and b edges 
covered by dimers in the thermodynamic limit. This leads to finding the expected number of 
creations or the expected number of annihilations. These provide an insight into the different 
behavior of the model depending on the value of u. Using the expectations computed in this 
section, we are able to determine the behavior when x — > 0. We also provide a duality result 
between high and low anisotropy. In Section [5J we concentrate on the behavior of the model 
when u < u c is fixed (and independent of x) and prove Theorem [2j Section [5] provides the 
proof of Theorem 0] when 7 = — 1 and we give the proof of Theorem [6] . The case when 
u = Ui is considered in Section [61 and there, we show the proof that the dimer model when 
u = Ui is equal in distribution to the noisy voter model. This leads to the convergence to the 
Continuum Noisy Voter Model. In Section [7J we prove the rest of Theorem 0] (for 7 7^ 1,-1) 
and also find an explicit expression for the correlation lengths for all values of 7. Section [8] 
gives some concluding remarks and directions for future work 

Acknowledgements: I am very grateful for the help of Richard Kenyon, who provided 
his valuable insight through countless hours of discussion. I would like to thank David 
Brydges for many discussions on statistical mechanical models, Jon Warren for pointing 
out the Brownian Web references, Zhongyang Li for comments on this paper, and Richard 
Arratia for sending a copy of his PhD thesis, [5rr79]. 

2. The Dimer Model on The Fisher Lattice 

In this section, we introduce some of the theory of the dimer models. Our brief survey 
contains only the most relevant results crucial for the rest of the paper. For a full treatment 
of results on dimer models on non-bipartite graphs see [BdT08| . We provide some analysis 
of the characteristic polynomial (defined in Section 12.21) . which sheds some light on the 
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behavior of the particle model. Finally, we give some facts regarding Pfaffians including an 
inclusion-exclusion formula. 

2.1. Thermodynamic Limit. Let G be a planar graph with structure given by the Fisher 
lattice, see Figure CD Suppose that G has edges weights x, ux or 1 depending on whether 
we have a vertical edge, a non-vertical edge off the triangle decoration or an edge on the 
triangle decoration. G is Z 2 -periodic meaning that the action of translations in Z 2 is an 
automorphism of G. Let G n be the quotient of G by the action raZ 2 . In other words, G n is 
a finite graph, with Fisher lattice structure, on a torus. 

A dimer is an edge and a dimer configuration, M, is a subset of edges of G such that each 
vertex is incident to exactly one edge. For a dimer configuration, we say that a 'dimer covers 
an edge' if that edge is present in the dimer configuration. Let Ai(G) denote the set of all 
dimer configurations of the graph G and let M(G n ) denote the set of all dimer configurations 
of the graph G n . On Ai(G n ), define the probability measure \x n for M G M(G n ) by 

. , (ux) Nb x Na 

where N a = N a (M) and Nb(M) are the numbers of dimers covering a and b edges in the 
dimer configuration M and Z(G n ) is the partition function: Z(G n ) = YlMeM(G n )^ ux ) a ' 
A measure of the form in (|2.ip is called a Boltzmann measure. 

A Gibbs measure on M{G) has the following property: Fixing the perfect matchings inside 
an annular region, then the perfect matchings inside and outside of this annular region are 
independent. Let T be the cr-field generated by cylinders (here, a cylinder is the set of dimer 
configurations of G containing a fixed finite subset of edges of G). Then, Theorem 6 of 
[LMT08] guarantees the existence of a unique probability Gibbs measure \x = fi(u, x) defined 
on (A4(G), J-) such that it is weak limit of the Boltzmann measures \i n as n — > oo for any 
dimer model on a non-bipartite graph. This limiting measure is called the thermodynamic 
limit. We denote [j,q = lim.x-»o fj,(u,x), F^O) to the be the probability of the event (•) under 
\i and E„ to be the expectation under \x. 

2.2. Partition Function, Local Probabilities and Notation. Here, we explain some of 
the theory of dimer models. In particular, we state results that are useful in the rest of the 
paper, such as, the formula for the partition function and local statistics formula. For proofs 
of these results, see |BdT08] . 

We refer to the graph G\ as the fundamental domain. Orient the edges of G\ so that the 
number of counter clockwise edges per face is odd. This is called the Kastelyn orientation. 
Such an orientation is given in Figure [6] for the fundamental domain. Assume that G (the 
infinite planar graph) has the periodic Kastelyn orientation induced by G\. As G is Z 2 
periodic, we can write any vertex of G as a vertex of Gi plus a translation of the fundamental 
domain, where the first co-ordinate of the translation is in the direction of the vector (1, 1) 
and the second co-ordinate is in the direction of the vector (—1,1) (i.e. v + (x,y) is the 
vertex v in the fundamental domain translated by (x — y, x + y)). 

Let V{G\) = {vi} be the set of vertices of G\ labelled as per Figure [U Let 71 (and resp. 
72) be a path in the dual of G%, winding around the torus in the direction of the vector 
(1,1) (and (-1,1)). Write e^- for the weight of the edge ViVj. For parameters z and w, let 
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Figure 6: Fundamental domain. The edge from v\ to vq crosses 71 and the edge from V2 to 
«5 crosses 72. 

Wij = z ril w n2 eij, where 

(1 Viv) crosses 71(72) from below 

ViVj does not cross 71(72) 

—1 V{Vj crosses 71(72) from above 

Let K(z, w) be the 6x6 anti-symmetric matrix given by 

Wij if Vi ~ Vj, and Vi — > Vj 



(2.3) 



(K(z,w)) 



—Wij if Vi ~ Vj, and Vj 
otherwise 



where Vi ~ Vj means that Vi and Vj share an edge and Vi — > Vj denotes an arrow from Vi to 
Vj in the Kastelyn orientation. This gives 



(2.4) 



K(z,w) 



I 
1 

-1 



y -bz 



-1 1 b/z \ 

-1 b/w 

1 a 
-a 1-1 

-bw 0-10 1 

1-10/ 



where a = x and b = ux. 
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The relevance of the above discussion is that we can write the partition function and 
probabilities of observing a local configuration in the thermodynamic limit as expressions 
involving K(z,w) (see [BdT08] and [KOS06J). The characteristic polynomial for the funda- 
mental domain, P(z,w), is defined to be the determinant of K(z,w). This gives 

(2.5) P(z, w) = a 2 + 2b 2 + a 2 b 4 + ab(l - b 2 ){z + - + w + -) + b 2 (l - a 2 ){- + -) 

z w w z 

The logarithm of the partition function can be written in terms of the characteristic poly- 
nomial, namely 

(2-6) l °Z Z = -^hs[ [ log^^) — - 

2(4tH) J ]z]=1 J\ w]=1 w z 

Let E = {ei = u\U2, . . . ,e m = U2m-iUm} be a subset of edges, with each iij representing a 
distinct vertex. }BdT08j showed that the probability of observing E in thermodynamic limit 
is given by 

(2.7) P M (ei, . . . ,e m ) = {\\ K ( U 2i: «2i+i)J Pf ({K~ l {iii, ttj)) T )i<ij<2m 

where T represents the transpose of a matrix and assuming v and v are in the same funda- 
mental domain, 

(2.8) K -\ v ^ + {x,y)) = -^ [ K-\z,w) v ^w x zy — -. 

This formula was discovered in |Ken97j for bipartite graphs. (|2.T[) is referred to as the local 
statistics formula. 

For the rest of the paper, we use the following notation. We denote Vi(n, m) = v % + (n, m) 
where i € {1, 6}. In other words, Vi(n, m) is the vertex Vi in the (n, m) fundamental domain. 
Note that u» = «i(0,0) = ^ + (0,0). 

2.3. Critical Values. In this subsection, we allow u £ (0,oo). A dimer model is said to be 
critical if there exists values of (z, w) S T 2 such that P(z, w) = 0. The dimer model on the 
Fisher lattice with weights (u, x) is critical when (—1 + vT-j- x 2 )/x 2 and (1 + vT+ x 2 )/x 2 . 
By Lemma [3 the particle model has the same distribution at each of these values of u. We 
define u c to be (—1 + vT+ x 2 )/x 2 , which is approximately 1/2 — x 2 /8. The dimer model on 
the Fisher lattice has other values of interest, namely (1 — y/l — x 2 )/x 2 and (1 + Vl — x 2 )/x 2 . 
Again, the particle model has the same distribution at each of these values of u by Lemma 
[71 We define Ui to be (1 — y/l — x 2 )/x 2 and refer to it as the independent value of u. The 
details of u = Uj and why it forms an 'independent' model is contained in Section [6l 

2.4. Particle Model probability space. In this subsection, we define the underlying prob- 
ability space for the particle model by use of a certain metric. We use the same space defined 
in |AB99| and [Cam08] . Let M 2 denote M 2 U {oo}, i.e. the compactification of M 2 . Let d be 
the spherical metric, which is defined to be d : K 2 x M 2 — > R by 

(2.9) d(u,v) = inf J(l + \<t)\ 2 y l ds 

where the infimum is over all smooth curves connecting u and v, <f) is parameterized by the 
arc length s and | • | is the Euclidean metric. 

Let H. be the space of all collections of curves in R 2 . We need to define a metric on % in 
order to determine its open sets. Let 71 and 72 be two curves with their parametrization at 
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time t given by 71 (t) and 72 (t) . Consider a complete separable metric space S of continuous 
curves in M?, with distance given by 

(2.10) 0(71,72) = mf sup d(7i(t), 72 (t)) 

te[o,i] 

where the infimum is over all choices of parametrizations of 71 and 72. The collection of all 
closed sets of S is exactly %. Let F\ and Ti denote two closed sets of curves. Then, D 
induces the Hausdorff metric, Dy, defined by 

(2.11) D^i{Fi,F2)) < e (V71 G ^1,372 G T% such that 0(71,72) < e and vice versa) 

Therefore, (%,-D%) defines a complete separable metric space. Let J-y_ denote the Borel a 
field associated with the metric Dy. Let My. be the space of probability measures taking 
values in It is clear that the particle model (before rescaling) are measures in the 

space M-u for all values of u and x. 

We remind the reader that we say the scaling window (x a ,x^) of some measure is the 
limiting measure obtained when taking the box (x a , x 13 ) to (1, 1) while simultaneously sending 
x — > 0. It is evident that the CNVM is in M% because the Brownian Web (through which 
the CNVM is contructed) is a random variable taking values in (T-L,J--u). For each u 7 and 
7 ^ 1, it is expected that the limiting measure of the particle model in the scaling window 
(x, x 2 ) is in M.%. 

For u < u c (u is independent of x), it is enough to consider the space .Mr2, which is the 
space of locally finite point processes on IR 2 , for the particle model. This is due to showing, 
in Section^ that there are no infinite lines and the loops formed by the trajectories of pairs 
of particles are finite P^-almost surely. 

For each n 7 and 7 / 1, the distribution of the locations of particles along an arbitrary 
horizontal line is contained in .Mir, the space of locally finite point process on M. By choosing 
the scaling window so that the density of particles along a horizontal line remains constant, 
the distribution of particles along this horizontal line converges to some measure in .Mr. 
This limiting measure, under the appropriate re-scaling, is found in Theorem [H 

2.5. Distribution of many particles. This self-contained subsection focuses on a formula 
for the locations of particles in the particle model. Essentially, this formula is an inclusion- 
exclusion type formula. We prove it using the recursive definition of the Pfafnan. However, 
another proof of the formula can be found using the fact that a Pfafnan of an anti-symmetric 
matrix can be thought of counting perfect matchings of a weighted graph (negative weights 
are allowed here), see [God93] for more details. 

The notation is set up to avoid explicitly writing out complete matrices. For 1 < % < n, 
let Xi represent the existence of a particle at Xj. This is equivalent to the appropriate a 
edge not being covered by a dimer in the underlying dimer configuration. The complement 
of Xi shall be denoted Xf, which represents a dimer covering the appropriate a edge in 
the underlying dimer configuration. Let P denote the law of observing the particles. Let 
K n (x\, . . . ,x n ) = {ki.j} 2r j =1 represent the 2n by 2n anti-symmetric matrix with fc2i-l,2i 
1 = W(Xf) for 1 < i < n and kij denote the appropriate inverse Kastelyn entries so that 

(2.12) P(X 1 C , ...,X c n ) = Pf (K n ( Xl , . . . , x n )) 

Define, K n [x\, . . . ,x k ,x k+1 , . . . ,x n ) = {&i,j} 2 ™ =1 , with k 2 i-\^i = — 1+v for 1 < i < k and 
kij = hj. Let I„ k = {aj.j} 2 ™ =1 represent the antisymmetric matrix with a,2i-i,2i = 1 for 
i £ {1, . . . , n) and aij = otherwise. Then 
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Lemma 3. 

(2.13) P{X 1 ,...,X n )=Pi(I n k -K n {x 1 ,...,x n )) = (-l) n Pi(K n { Xl ,...,x n )) 

The proof does not require the edges covered to have the same probability of being covered. 
The result is stated in this fashion for convenience later in the paper. 

Proof. The proof follows by using an iteration, with the key step relying on using the recursive 
definition of the Pfaffian. Consider 

(2.14) P(X U XI ...,X c n )= F(Xl ...,X c n )- ¥(Xt, . . . , X n ) 

with P(X|, ...,X%) = Pf(iT n _i(x2, . . . , x n )) and P(Xf, ...,X°) is given in (j2T2]l . On the 
other hand, by using the recursive definition of the Pfaffian and. letting A — K n {x\^ X21 • • • •> ^n) 

2n 

(2.15) Pf {K n (x u x 2 , • • • , x n )) = (-1 + v)Pi(A iA ) + J^-l^Pf {A^ 

i=3 

where A; j is the matrix obtained by removing both the 1st and i th row and column. Notice 

that Pf(A M ) is exactly Pf {K n ^(x 2 , . . . , x n )) while uPf(A ij2 ) + £-= 3 (- 1)^1,^(^1,;) = 
Pf(K n (x\,X2, ■ ■ ■ ,x n )). Therefore, substituting (|2.15|) into f)2 . 14|) gives 

(2.16) ¥(X U XI ...,X c n ) = -Pi(K n ( Xl ,x 2 , • • • , x n ))- 
Applying the steps to obtain (|2.16p iteratively gives the result. 

□ 



3. Thermodynamic Limit 

This section focuses on the behavior of the model in the thermodynamic limit. We in- 
troduce a generic method for calculating expectations, which is shown by computing the 
expected number of dimers covering b edges per fundamental domain. This method can be 
used to compute other expectations, such as the expected number of dimers covering a edges 
per fundamental domain. We then use expansions of the computed expectations to find the 
behavior of the model in the limit x — > 0. All the results are reliant on analysis of P(z,w), 
which we provide later in the section. Finally, we state and prove a 'duality' statement for 
low and high anisotropy. 



3.1. Expectations for < x < 1 and < u < 1. Let Nj, denote the number of dimers 
covering b edges per fundamental domain. By the setup of the model, there are 2 b edges in 
each fundamental domain. Let N a c denote the number of a edges not covered by a dimer per 
fundamental domain. Let Nx denote the number of creation singularities per fundamental 
domain. In this subsection, we find exact results for E^Ay, E,^[N a c] and E^[iVx]- Figure [7] 
shows a plot of E[iV~b] and E[JV ] for x = 0.5. Recall that we defined u c = (— 1 + yl + x 2 )/x 2 . 



Theorem 5. For < x < 1 and < u < 1 



(3.1) 



/(*,«, i) + 

2v?x 2 



ifu<u c 
ifu = u c 

k f(x,u,l/r) + JtXSx* i f u>u c 
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Figure 7: A plot of E[A^] and E[iV a ] against u G (0, 1) for x = 0.5. Notice that the difference 
between the two curves is 2E[A r x]- 



where f(x,u,r) is a continuous function given by 
(3.2) 

f{x, u, r) = C{x, u)((u 2 (-l + x 2 ) + rV(-l + x 2 ) + r(l + uV - 2u 2 (l + 2x 2 )))/C 



rs — 1 



(3.3) + c(x )U )(-i + rV(-i + ^)(n('^,^4V(n^ s r s 



rs — 1 rs — 1 / \rs-l rs-1 



/or r ^ 1 where 

(3.4) v^(l + u 2 * 2 ) 

(3.5) r = 



irr(rs — l)u 2 {— 1 + x)(l + a?)(— 1 + it 2 x 2 ) ' 
1 + 2u 2 - u A x 4 + (-1 + uVyi - 4n 2 + 2u 2 x 2 + u 4 x 4 



2(— it — 2u 2 x — u 2 x 2 ) 



(3.6) 



1 - +2u 2 - u A x A + (-1 + n 2 x 2 )Vl - 4u 2 + 2u 2 x 2 + u 4 x 4 
2(— it + 2u 2 x — u 2 x 2 ) 



fC{k) represents the complete elliptic integral of the first kind with modulus k and U(n,k) 
represents the complete elliptic of the third kind with characteristic n and modulus k. Also, 

(-1 - x 2 + yT+gK-2 - x 2 + 2yTT?) / 1 x + 

(3.7 /(x,it,l) = — -—log r + l 

(-1 + VI + x 2 ) 3 \ 2iri x-i 

The elliptic integrals are defined for u > u c by an analytic continuation argument (see 
[Law89]). Note that E^fiVfc] is a continuous function. 

Proof. There are two ways to compute E^[iV&]. Either use the inverse Kasteleyn matrix to 
compute the probability of seeing an edge along u 5 to t>2 (and multiply the corresponding 
result by 2) or we can use the properties of the partition function. We follow the latter 
approach. Differentiating under the integral sign in (|2.6p with respect to b and multiplying 



16 



SUNIL CHHITA 



by b gives 

(3.8) MJVJ-A/ / A^r-- 

8vr 2 y| w | =1 i| 2 | = i P(z,u>) z tu 
Taking the change of variable z = vw gives 

(3.9) W1 = A/ / i^)*^ 

8tt 2 y [w | =1 J| Z | =1 P(vw,w) v w 

Setting a = x, b = ux, then the roots of the denominator of the integrand with respect to 
w are given by 0, r\(v) and r2(u), where 

(3.10) 

1 2 2 22,22,222 4 4 At \\ 

v — 2u v — u v + u x + u v x — u vx =F A[v) ) 



-u 



2(uv + uv 2 — u 3 vx 2 — U 3 V 2 X 2 ) 
where 

(3.11) A 2 (v) = (-4(n + uv - u 3 x 2 - u 3 vx 2 )(uv + uv 2 - u 3 vx 2 - u 3 v 2 x 2 )+ 

( 2 , , o 2 , 2 2 22 222,4 4\2\\ 

[u +v + 2uv + uv — U X — U V X + u vx ) )). 
The root r\(v), lies within the region \w\ < 1 for all \v\ < 1 while ^(v), lies in the region 



\w\ > 1 for all \v\ < 1. Therefore, the first integral in (|3.9p can be evaluated using residue 
calculus. 

The residue at w = for (|3,9p is given by 



ux 1 /" -l + 3u 2 x 2 , l-2>u 2 x 2 
( 3 - 12 ) ~o~ oZ7 / —< i , TSZST ^ 



2 27rz 7|„i = i uvx(— 1 + u 2 x 2 ) 2(1 — tt 2 x 2 ) 
The residue at w = ri(f) for (|3.9p is more complicated and is given by 

27ri7 w=1 2?;(-l + n 2 x 2 )^(T;) 

where 

(3.14) Q(v) = (1 + uV)(v + u 4 ra 4 + u 2 (-l + x 2 + w 2 (-l + x 2 ) - 2u(l + 2x 2 ))) 
The residue at v = for f)3. 13|) is given by 

1 + u 2 x 2 

(3-1^) ; 2 — 2T" 

—2(1 + u z x z ) 

Notice that the contribution from f)3. 12|) and f)3. 15|) is G(x 2 ) and gives the additive term 
shown in the formula for Eu[iVjJ. It remains to finish evaluating (|3.13|) . which is done by 
analysis of the roots of A(v). The roots of A(v) are r, s, 1/r, 1/s, where r and s are defined 
in Theorem [3 The behavior of the underlying dimer model is dependent on the location of 
these four roots, which is given explicitly in Section 13.31 

Notice that \s\ > 1 for all values of u. However, |r| > 1 if u < u c , r = 1 if u = u c and 
\r\ < 1 if it > u c . We can consider each case separately. Note that when u = Ui > u c , 
A(v) breaks down into a quadratic. In this special case, we can calculate using residue 
calculations. 

For u < u c , the contour of integration in (|3.13p can be deformed to a contour integral 
surrounding the branch cut from 1/r to 1/s as the residue at v = has already been 
considered. This integral can then be deformed further giving two line integrals, which can 
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both be written as integrals from 1/r to 1/s. These integrals are given by elliptic integrals, 
which can be simplified to the desired form. 
For u = u c , then A(v ) is given by 



(3.16) (v - l)y/{v - s)(v - 1/s). 

This implies that the contour of integration in (|3.13p can be deformed to a contour integral 
surrounding 1 and 1/s. Completing the integrals via line integrals gives the required result. 

For u > u c , then we can repeat the same steps as given for u < u c but with the branch 
cut running from r and 1/s. 

□ 

Corollary 1. For < x < 1 and < u < 1, 

2 

g{x,u,r)-f^z ifu<u c 



(3.17) E^[N a 



2 

g(x, u, 1) - ifu = u c 

2 

9{x, u, 1 /r) - j^j ifu>u c 



where 
(3.18) 

g(x, u, r) = D(x, u)(u 2 (l + x 2 ) + rV(l + x 2 ) - r(l - 2u 2 + u A x 4 ))K 



(r 



(-1 + rs) 



(3.19) i% u )(-i + rV(i + ^)(n(^,^-(n^ s r 



rs 



1 ' rs — 1 / \ rs — 1 ' rs — 1 



forr^ 1, where D(x,u) = C{x, u)(— l+u 2 x 2 )/(l+u 2 x 2 ), C(x,u), r, s, K andU are defined 
in Theorem^ Finally, 

(3-20) g(x, u, 1) = ( 1 " 1 lo S (— 

1 — x z \ m \x — 

Proof. E,^[N a c] can either be computed using the partition function or by using inverse 
Kastelyn techniques, namely 

(3.21) E M [iV>] = l-xK~\v 3 ,v±). 

This can be computed using the same techniques used in Theorem □ 

A creation consists of 3 edges, so its expectation consists of a Pfaffian of a 12 by 12 matrix. 
This is somewhat messy. However, the next lemma shows that E^JiVx] can be computed 
using E^Nb] and E^iVj. 

Lemma 4. 

(3.22) E^Nx] = I (E^Nt] - E^iV]) 

Proof. The trajectories of paths in the particle model form loops and infinite lines. These can 
be represented on the underlying dimer model as dimers covering b edges and not covering 
a edges. Let Nb denote the number of dimers covering b edges, which do not belong to 
creations, per fundamental domain. Clearly, for each loop or line, Nb = N a c and 

(3.23) EJN b )=EJN a c]. 
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On the other hand, 2Nx = Nj, — iVj, because for each loop or infinite line, the difference 
Nb — Nfj is the number of pairs of b edges seen in a fundamental domain (per fundamental 
domain). These are exactly creations (from definition of Nf,). Therefore, 



(3.24) E^[N b ]-E^[N b ] ■. 

Combining (|3.23[) and (|3.24p gives the result. 

This leads to the following 
Corollary 2. For < x < 1 and < u < 1 

f(x,u,r) - g{x,u,r) 



2EJN 



□ 



(3.25) 



EjiV 



x 



f(x,u,l) 
f(x,u,r) 



g(x,u, 1) 
g(x,u,r) 



x 2 (l-2u 2 +u 2 x 2 ) ' 

(1-X 2 )(1-U 2 X 2 ) 

x 2 (l-2u 2 +u 2 x 2 ) '' 
(l-x 2 )(l-u 2 x 2 ) 
x 2 (l-2u 2 +u 2 x 2 y 
(l~x 2 )(l-u 2 x 2 ) 



if u < u c 
if u = u c 
if u > u c 



where the functions f and g are defined in Theorem \3.9\ and Corollary [0 

3.2. Expectations in the limit x — > 0. The results in the previous subsection can be 
expanded out in terms of x for some values of u. By computing the expansions of E^fiVjJ 
and E/JiVx], we can determine the behavior of the model in the limit x — >• 0. This does 
require some attention - taking the limit through the underlying integrals cannot always be 
done. Therefore, we need careful expansion of the appropriate elliptic integrals. 

Lemma 5. 

0(x 2 ) if u < u c (fixed) 

f + 0(x 2 ) ifu = u c ' 

x + 0(x 2 ) if u = m 

2 - 2 1 + 0(x) ifu>Ui (fixed) 

where 6 = arccos(— l/(2u)). Furthermore, the first order term o/E At [A r a c] is equal to the first 
order term ofE^Nb] in each case. 

Note that 2/tt < E[N ac ] < x and E^[JV a c] is increasing for u c < u < u\ by (|7.5p . The same 
is true for E^fiVj]. 

Proof. For u < u c (fixed), the modulus of the elliptic integrals, k = (s — r)/{rs — 1) = 
4it 2 x/\/l — 4u 2 + 0(x). The expansions of the elliptic integrals are given by 



(3.26) 



(3.27) 
and 
(3.28) 



U(rk,k) -U(k/r,k) 



fc7r(-l + r 2 



+ 



3/t 2 vr(-l + r 4 



+ 0(k 3 



4r 16r 2 
Expanding r and s in terms of x and simplifying gives the result. 

The expansion for u = u c follows by a series expansion of the expression given in Theorem 

m 

For u = Ui, we can calculate E„[iV&] using solely reside calculations. This is due to the 
branch cuts contract to a point. Hence, A{v) defined in the proof of Theorem[5]is a quadratic. 
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For u > Ui (fixed), the modulus of each of the elliptic integrals is given by i/k' where 
k' = \s/r — l\/\s — l/r\ = Au 2 x/\/Au 2 — 1 + 0(x) and i = yf—1. The expansions for the 
elliptic integrals are given by 

(3.29) K{k) = k! log (y J + k' log 2 + 0(k' 2 ) 

and 



ir i \ „ / i i \ 1 



„'2\ 



(3-30) n {T-J?)- n {^v)= k ' l °Z7 + °( k 

which are both computed in the Appendix of [Chhllj . Plugging back the relevant expressions 
gives the result. 

The last statement follows by finding the expansions E^[iV a c] which is achieved using the 
same method. 

□ 

The expansion for E^fiVx] can be also computed for certain values of u. This is given by 
the following 



Lemma 6. 



(3.31) E^lNx] 



1(1- 2u 2 - Vl - 4u 2 )x 2 + 0(x 3 ) ifu<u c (fixed) 



2 y 

^+0(x 3 ) ifu = u c 

\ + Q(x 3 ) ifu = Ui 

(— ^\/4u 2 — l)a; 2 logx + 0(x 2 ) if u > Ui (fixed) 



Proof. For u = or u = u c , the exact same method employed in the proof of Lemma [5] 
holds. Otherwise, notice that E^JiVx] is of the form 

(3.32) Cl (x,u)K(k) + c 2 (x,u)(U(r/k,k) - IL(l/(rk), k) + 0(x 2 ) 

where c\(x, u),C2(x, u) represent coefficients of the complete elliptic integrals and k represents 
the appropriate modulus of the elliptic integrals, given in Theorem [5l For u < u c , c\ and C2 
are both of 0(x 2 ), hence the expansions of the elliptic integrals given in the proof of Lemma 
[5] hold for u < u c . When u > Ui, ci(x,u) are both O(x), which implies that we can use the 
expansions for the elliptic integrals in their given form, without needing to compute extra 
terms. □ 

We can now prove Theorem [TJ 

Proof of TheoremUl For u £ (0,1), using Markov's inequality and the Borel-Cantelli Lem- 
mas, we obtain 

(3.33) lim FJN X = 0) = 1. 

x— >0 

As the distribution of the creations is the same as the distribution of the annihilations, (|3.33p 
implies with probability 1, there are no annihilations or creations for u G (0, 1). By Lemma 
in the limit x — > and u > 1/2, particles have density 2 — 20/tt. For u < 1/2, by applying 
Markov's inequality and the Borel-Cantelli Lemmas, 



(3.34) hm P M (iV ac = 0) 



□ 
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3.3. Analysis of P(z,w). This subsection is primarily a consequence of the calculation in 
Theorem^ We provide the behavior of the roots of A(y), defined in (|3.1ip . As a consequence, 
we can prove a 'duality' result, which gives an isometry for high and low anisotropy. 

Recall that P(z,w) = det K (z, w) represents the characteristic polynomial and that a 
dimer model is said to be critical if P(z, w) = for some (z, w) S T 2 . Recall that r, s, 1 /r, 1 / s 
are the roots of A(v), as defined in (|3.1ip . The values of the roots of A(v) are dependent on 
the anisotropy as follows: 

• For u < u c , the roots are r, s, 1/r and 1/s for r,sfR with 1 < r < s. 

• For u = u c , the roots are s, 1 (twice) and 1/s for sGM with 1 < s. 

• For u c < u < Ui, the roots are r, s, 1/r and 1/s for r, s £ R with r < 1 < s. 

• For u = Ui, the roots are r (twice) and 1/r (twice) for r£K with r < 1. 

• For u > Ui (u < 1/x), the roots are r,r, 1/r and 1/r for r £ C with |r| < 1. 

It is evident that for u = u c , A{— 1) = 0, which implies the model is critical when u = u c . 
Taking u = Ui, then all entries in the inverse Kastelyn matrix can be computed solely using 
residue calculations. This leads to the model possessing independent trajectories of particles 
(see Section [6] for further details) . 

Provided that u £ (0, 1), at low temperature the distance between certain pairs of roots 
is O(x). For u < u%, r — s is 0(x). For u > Ui, r — 1/r is 0{x). Notice that choosing 
u = 1/2 + x a c for < a < 2 and c is a constant, r — 1/r is 0{x a / 2 ). 

The next lemma gives a 'duality' for ux. Here, we temporarily drop the assumption that 
u £ (0, 1). Let X denote any configuration of a edges not covered by dimers in the dimer 
configuration. For notational simplicity, write P M ( UiX ) = f(u,x,X). 

Lemma 7. For u < 1/x and x > 0, 

(3.35) F(u,x,X)='p(—f,x,X) 

This shows that the particle model has the same distribution when choosing the anisotropy 
to be either l/(ux 2 ) or u, provided u < 1/x. This also implies that the two critical points 
(— 1 + vT+ x 2 )/x 2 and (1 + y/l + x 2 )/x 2 are equivalent (in distribution) in terms of the 
particle model as (1 + vT+ x 2 )/x 2 = l/(u c x 2 ). 

Proof. Let K~ 1 (vi,Vj(n,m))(u, x) denote K~ 1 (vi,Vj(n,m)) with parameters u and x, where 
i,j £ {3, 4}. By following the techniques in the proof of Theorem [5] and only computing the 
integral with respect to w gives 

(3.36) ^ 1 (f i ,u j (n,m))(l/(nx 2 ),a;) = (-l) n+m K~ x Vj (n,m))(u,x). 

In other words, we have an two integrals with respect to v which differ by a factor of (— l) n+m . 
Taking the Pfaffian of the matrix whose entries are given by the required inverse Kastelyn 
entries gives the required result (the factor (— l) ra+m always affects an even number of rows 
and columns). □ 

4. u < u c (not dependent on x): 

In this section, we prove Theorem [2] using a sequence of lemmas and propositions. These 
lemmas and propositions lead to the following informal description of the model. Before 
re-scaling, particles are created in pairs and the probability of survival of these particles 
decays exponentially with time. It is not unreasonable to expect that the particle model 
only consists of loops. The interactions between each of these loops can be shown decay 



PARTICLE SYSTEMS ARISING FROM AN ANTI-FERROMAGNETIC ISING MODEL 21 

exponentially with distance and with time. Under any scaling window (x, x) each loop 
contracts to a point and the distribution of these points is a Poisson point process with 
intensity v?{l/yj\ — 4u 2 — 1) as given in Theorem [2l Figure [8] shows a realization of the 
particle model for u < u c . 



Figure 8: A simulation of the particle model on a grid of 100 by 100, with u = 0.4 and 
x = 0.1 using Glauber dynamics 

Before embarking on stating and proving the lemmas and propositions, we first introduce 
and recall some notation. Recall that Vi(m, n), . . . , VQ(m, n) represent the vertices in the 
fundamental domain (m, n) and that the fundamental domain (m, n) is the fundamental 
domain at (0,0) translated by (m + n,m — n) G Z 2 . Define X(m,n) to be the event of 
observing the leftward upward b dimer leaving the fundamental domain (m, n) (i.e. the 
event of seeing a dimer covering the edge from Ve(m, n) to v\(m + l,n)). The inverse 
Kastelyn entries for vertices in different fundamental domains requires Fourier coefficients 
of P(z,w) denoted by H(m,n). In other words, 

Iff z m w n dz dw 

(4.1) H(m,n] 



(2ttz) 2 _/ n=1 _/ h=1 P( z ,w) z w ' 
The next lemma gives H(m, n) as an expansion in terms of x. 

Lemma 8. Let T\ and be the roots of the polynomial u + w + uw 2 with \r\\ < 1 and 

1 7^2 1 > 1 and set s% = u/(u + r{) and S2 = u/(u + r^). Then, for n,m > 1, 

(4-2) 

(- U )m+n (m + n)! ^(1,-m.l + n,^) ^ {-s 2 ) n - 1 - k \ n( 

H(m, n) = r — - S 2 > 777 7-tt + 0(x) 

u (r l — r 2 j \ m\n\ kl{m + n — ky. J 

where 2-F1 is the confluent hyper geometric function of the second kind. 
Proof. By letting 

(4.3) P(z, w) = lim — zw det K (z, w) = (1 + zu + wu)(wz + uz + uw) 
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then using the analysis of the roots in Section 13.31 by the Dominated Convergence theorem, 
we can take series expansion in terms of x, which yields 

Iff z m w n 
(4.4) H(m,n) = -—2 / w -dzdw + 0(x). 

( 27 "r J\z\=iJ\ w \=i P(z,w) 

As |ttu;|/|ii + w\ < 1 for all \w\ = 1, then using the factorization of P(z,w) in fj4.3|) gives 

1 f ( uw 
2m J\ w \=i V u + w J u + w + uw' 2 



(4.5) H(m,n) = — / ( -— — ) - ; - ; 2 dw. 



Notice that 



2 



1 1 



(4.6) u + w + uw , 

U{r\ — T2) \W — n W — T2 

As |ri| < 1, the contribution for (|4.5p obtained by deforming the contour of integration to 
\w — r\\ = e is given by 

(_ v )n m+n 

(4.7) 1 U} ri 



u{r2 — r\){r\ + u) n 

The remaining contribution for (|4.5j) comes from deforming the contour of integration to 
\w + u\ = e. For ti£C with v ^ u, notice that 

(4.8) 

i I /_^y^ dm m+ g /_u_y-» 

27T« 7|tu+u|=e V U + WJ W + V kl{m + U - k)\ \U - V J 

(4.9) = -(-l) m f — ^ v m+n 

\ u — V J 

(m + n)! 

m!n 

Substituting in the desired values of v and manipulating gives the desired result. 



(4.10) - {-u) m+n y 2-Fi(l, -m, 1 + n, 1 - 



□ 



From the above lemma, we can compute the asymptotic expansion of H(m, n) for m, n > 0. 
Let m = n + k, then 

(4.11) H(n + k,n)= (-^(^ + ^' + 0(t| »^f 1) 

We also require H(—m,n) or H(m, —n) for m, n > 1, which is given below. 
Lemma 9. For m,n> 1, then 

(, \ n.+m 
^ +0(X) - 

Proof. Without loss of generality, suppose that n > m. Taking the series expansion about 
x = gives 

Iff w n 

(4.13) H(-m,n) = - ^ / / ; -dzdw + 0{x) 

(2my J\ w \=i J\z\=i z m {l + zu + wu)[uw + uz + wz) 
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By noticing that 

(4-14) 7 ^ r = J U + W 

(1 + zu + wu)(uw + UZ + WZ) U + W + UW 1 \UW + UZ + wz 1 + uw + uz 

then completing the residue calculations with respect to z leads to the simplification 

1 r (-l\ m v m w n 

(4.15) H(-m,n) = / i , —^- 7 — dw + 0(x) 

v 1 y ' ' 2m J\ w \ =l (u + w + uw 2 ){l + uw) m v; 

(4.16) = { } V 1 +Q(x) 

(1 + uri) m Vl - 4u 2 

where ri is the root lying inside of the unit circle of the polynomial u + w + tiro 2 . Simplifying 
the above equation gives the result. □ 

It is clear that for u < 1/2, then H(m, n) decays exponentially. This leads to the following 
proposition: 



Propostion 1. For m,n 6 Z and \m\ > \n\ with k = \ \m\ — \n\ 

\A-4u 2 , 



(1.17) F(X(0,0),X(m,n)) = 4uV ( 1 - -=1=^ + Cx 2 (-u) in+2k + 0(u 4n+2/c+1 ) 



where C is an arbitrary constant. 

Proof. We only compute the result for m > n > 1 as the other cases follow from very similar 
calculations. Using the local statistics formula (|2.7|) . P(X(0, Q),X(m, n)) is given by 

(4.18) 

/0 K^ 1 (v e ,v 1 (—1,0)) K-^v^VQi-m,-^) K~\vz, «i(-m + 1, -»)) \ 

2 2 pf ••• K _1 (wi,?j 6 (-m - 1, -n)) ■K" -1 (vi,vi(-m,-n)) 

U *T->6,«i(-l,0)) 

V ■■■ / 

Two of the interaction terms can be computed directly using Lemma [8] by noticing that 

(4.19) -ftT -1 (i>6, vq(— m, —n)) = —K~ l {v\,vi{—m, —n)) = —uH(m — 1, n) + uH(m + 1, n) 
It still remains to find K~ l (v\, vq(— m — 1, —ri)) and i^ _1 (f6, i> i(— m + 1, — n)). Notice that 

Iff z m ~ 1 w n 

(4.20) ^-1(^,^(^ + 1,-^) = ——-/ / ■ ■ dwdz + 0(l) 

x(2my J\ z \ =1 J\ w \=i uw + uz + wz 

< 42i » = i^;:?). '-'"' + "" +o ">- 

The remaining interaction entry, K~ l (v\, vq(— m — 1, —nj), contains terms of the form 

x f f z n w m 

(4.22) I(m,n) = -. -tt / / -. -7 -^dwdz 

(2lTl) z J\ z \ = l J\ w \ = i {uw + uz + wz){± + uz + uw) z 

I(m,n) can be computed by differentiating H{m,n) by n and ri- However, this leaves 
a contribution of 0(u 2n+k+1 ). Combining all the linear terms of I(m,n) and multiplying 
by K~ l (v§, vq(— m, — n)) contributes a maximum of 0(tt 4n+2fc+1 ), and so is absorbed into 
0(u An+2k+1 ). 
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The term K~ x (vq, «i(1, 1)) can computed via residue calculations, which gives 

(4.23) ux ( 1 ; 1 

V ' V VI -4n 2 , 

Putting together the various contributions gives the result. □ 

Proposition [T] can be generalized to show that the difference between the joint probability 
for k edges and the product of the probabilities of each edge is an exponentially decaying 
term. This can be achieved by using a complete graph representation of the Pfafhan and 
ignoring the matchings which contain exponentially decaying weighted edges. For brevity, 
we omit the details. 

Let A n be the event of seeing a string of b edges crossing an n by n square. This 
corresponds to seeing a particle traversing an n by n square in the particle model. Then 

Propostion 2. 

(4.24) P(,4 n ) < Ce~ Bn 
where B and C are positive constants. 

Proof. The event A n is contained in the event of seeing two b edges at the top and bottom 
of the n by n square at anys position. Therefore 

(4.25) F(A n ) < n 2 sup P(X(0, 0) n X(n - k, k)) 

k 

By a symmetrical consideration, the maximal probability occurs when the two b edges are 
vertical, i.e. when k = n/2. Using Proposition [T] gives 

(4.26) ¥(A n ) < n 2 (u 4 x 4 (l - - 1 ) + Cx 2 {-u) 2n 



Vl -4n 2 , 

Taking n ~ — logx shows that ¥(A n ) decays exponentially with distance. □ 

The above lemmas and propositions enables the proof of Theorem [2J 

Proof of Theorem By Proposition [21 each loop generated by a string of b edges is finite 
and does not appear in the scaling window by (x, x) almost surely. Therefore, each loop of b 
edges will be contracted to a single point under re-scaling. Tightness is guaranteed for point 
processes (see [DVJ88] ). It is enough to consider the locations of b edges and show that 
the finite dimensional distributions of the corresponding measures converge in the scaling 
window (x,x). We shall follow the approach developed in the proof of Theorem U] for the 
case u = u c (see Section 15. 2p and use the corresponding multi-index notation. 

Let ¥ x denote the point process in M 2 arising from the location of a single b edge in a 
loop. Let A4, for i £ {1, . . . , fc}, be disjoint, non-empty simply connected domains and let 
Ai be the union of lattice points in ((2Z x 2Z) U ((2Z + 1) x ((2Z + 1)) so that xA\ -> A { as 
x — > 0. For n € Nq, then 



(4.27) E x 



J2 ••• E nx(h,ti),...x(i w ,t ]n] )). 

(il,ti)EAi (i\ n \,t\„\)eA k 



As each edge is separated by a distance of at least then each of the interaction terms in 
the inverse Kastelyn matrix converges to zero due to the exponential convergence. Hence, 
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taking the above limit ELS X — > gives 
(4.28) 



/ . . . / u 2 H 



1 



Vl - 4u 2 - 1 



dx™ 1 . . . dx\ 



Let 
(4.29) 

and set 
(4.30) 



Qx(z) 



n 



N(Ai)\ 



\(N(A 3 )-p,)\ 



2|n| 



1 



Vl - 4u 2 - 1 



|n| 



i.e. a generating function for a Poisson point process. Notice that each generating function 
can also be written as a power series whose coefficients are given by the finite dimensional 
distribution probabilities for the corresponding measure on the sets Ai, . . . , A m . Following 
the same convergence argument given in the proof of Theorem [5] for u = u c , then it is 
clear that the derivatives of Q x (z) (with respect to z) converge to Qq{z). This proves the 
convergence of finite dimensional distributions. □ 



5. The particle model for u 



This section focuses on the critical case in which u = u c = - 
model is said to be critical if P(z, w) = for some (z, w) G T 2 
that the values of (z,w), which gives P(z,w 



Recall that a dimer 



Moreover, it can be shown 
0, must be real (see |BdT08^ . We first look 
at the local dynamics in the thermodynamic limit. This allows us to produce random walk 
type estimates that are required later. We can then prove Theorem[5]for u = u c (or 7 = —1) 
by obtaining the appropriate correlation kernel in the scaling window (x,x 2 ). Finally, using 
the random walk estimates, we can determine that the set of creations in the scaling window 
leading to 'non-trivial paths' (paths of length greater than e > in the scaling window) is 
a locally finite Poisson point process. Figure [9] shows a simulation of the particle model for 
x = 0.1 and u = u c . 

Our constraints for the choice of scaling window is that it satisfies Donsker's scaling (i.e. 
scaling space by the square root of time) and gives non-trivial paths. We choose to take the 
scaling window (x,x 2 ) by the following argument: Suppose we choose the scaling window 
(x a ,x 2a ). It suffices to look at the expected number of dimers covering b edges in a box of 
size (x~ a , x~ 2a ). By Lemma El this is G(x 1_3a ). The only value of a that ensures that we 
see non-trivial paths (i.e. expected number of dimers covering b edges is x~ 2a in the box of 



size (x 



-2a 



)) is a = 1. 



5.1. Local Dynamics. This section concentrates on computing random walk estimates of 
the particles in the particle model for u = u c . These estimates rely on computing the entries 
of an appropriate inverse Kastelyn matrix. The following lemma provides the underlying 
entries of the inverse Kastelyn entries. 



Lemma 10. For m > n and m = o(l/x) 
(5.1) 

1 f f z m w~ n - (-l) m - n dz dw 



(2vri) 2 



>|=l •/[*[=! 



P(z,w) 



z w 



Fo(m, n)+F±(m, n)+F2(m, n)+F^(m, n)+0(x 4 
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Figure 9: A simulation of the particle model on a 100 by 100 grid with u = u c and x = 0.1 
made using Glauber dynamics. 

where 



(5.2) F (m,n) = (-l) m ~"(m + n), 



, . . , 2(-l) m ~ n (m + n) 2 x 

(5.3) F 1 (m,n) = -± — - - — ; 

7T 

(5.4) F 2 (m, n) = --(-l) m " n (m + n) (l + 2m 2 + 4mn + 2n 2 ) x 2 

6 

and 
(5.5) 

/ i\m— n 

F 3 (m, n) = - (m+n)(-7-20m 2 +12m 4 +200mn+48m 3 n-20n 2 +72m 2 n 2 +48mn 3 +12n 4 

The expansion does hold for n > m. The (— l) m_n is required so that the above integral 
converges (due to having a singularity at z = w = —1 as P(— 1, — 1) = 0). The lower order 
terms are required for later computations involving the random walk estimates. 

Proof. The first integral can be computed via a residue calculation, leaving 



X 4 

(5.6) — f — 

2vTi 7^1 = 1 



[w|=l 2 2 + x 2 - 2V1 + x 2 K/(l + u>) 2 (1 + u> 2 + w (2 + 4x 2 )) 



Note that there is a branch cut situated along the negative real axis from 1 — 2x 2 — 
+ x 2 to r = —1 — 2x 2 + 2x\/l + x 2 . Therefore, the integral over |u;| = 1 in (|5.6[) can 
be deformed to an integral from w = —1 to w = r (above the branch cut), an integral from 
w = r to w = — 1 (below the branch cut) and a contour integral around w = \{x 2 — 1). 

The contributions given near the branch cut can be computed in a very similar fashion 
as follows. By taking a series expansion of the numerator about the point w = —1 gives an 
integral of the following form 
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(5.8) h := 1/(2™) I „ {W + 1 J 1 =dw 



(5.7) ± I* X \ f>- ^ + 

V ; 2m J_ x 2(2 + x 2 - 2VT+^) V( w + + w 2 + w(2 + Ax 2 )) 

where is a function of m, n and x. Each of the integrals 

r (w + 1)^ 1 

-1 y/(w + l) 2 + 4wx 2 

can be evaluated. It suffices to only compute Ij up to i = 5 as F = 0(x l ~ 1 ) and each is 
0(1). Summing up the contribution near the branch cut gives 

—Fo(m, n) + Fi(m, n) — F2(m, n)x 2 + F 3 {m, n)x 3 . 

It remains to calculate the other contribution in (15. 6j) , namely the contour integral around 
u> = \{x 2 — 1). Taking a series expansion in x gives 
(5.9) 

1 f 2(-l) m w m - n 2((-l) m w 1+m ~ n (l + 2w + m(l + w) 2 ))x 2 . 

' " -dw + 0{x ). 



2iri J\ w+ l\ =e (1 + 2w) m (l + w) 2 (1 + 2w) m+l (l + w) 4 

This can be computed using the following integral 

1 f w m-n 

2tu y| TO +i|= e (1 + 2w) m (l + w)« 

These can all be found by a higher dimensional residue formula and are given by G(l) = 
(-l) n ,G(2) = (-l) n (m + n) 

(5.11) G(3) = ^(-l) n (-2 + m 2 + n + n 2 + m(3 + 2n)) 
and 

(5.12) G(4) = -(-l) n (2 + m + n)(-6 + m 2 + n + n 2 + m(7 + 2n)). 

6 

Adding up the various contributions gives the result. □ 

Recall that Vi(n,m) denote the vertex Vi in the (n,m) fundamental domain. For (y,t) € 
(ZxZ)U((2Z+l) x (2Z+1)), let X(y, t) denote the event that there is a particle at (y, t) where 
we think of y as space and t as time. Note that ¥(X(y,t)) = 1 — P('a' edge) = 1 — xK~ 1 (3, 4). 

Propostion 3. The covariance of two particles a distance n away with n = o(l/x), is given 
by 

(5.13) F(X(0, t),X(n, t)) = ( 1 - ^] x 2 - —x 3 + 0(x 4 ) 

Proof. The proof of this proposition is a calculation using the inclusion exclusion formula 
given in LemmaH The interaction entries, K~ l (v$, vs(n, — n)) and K~ l {v 3 , v^(n, — n)), can 
be calculated as an expansion using Lemma [TOl 

te-xAS zr-l < / \ \ 1 8nx , , , 2 , 2 4(-3 + 34n + 32n 3 )x 3 

(5.14) K n t_ n (v 3 ,v 3 (n, -n)) = 1 + (1 + 4n 2 )x 2 - -i '— 

and 

1 / / 2 (10 + 8n)a; 2 o 

(5.15) X~i n (u 3) u 4 (n, -n)) = - - 2nx + - — - 3nx 3 . 



Pf 



x 2 /4 _^! + 0(x4) 

7T 
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In Section El we computed 1 — xK~ l (v 3 ,v 4 ). Plugging the required entries into the appro- 
priate Pfaffian gives the required result. □ 

The local dynamics of the particles can also be computed. This entails the conditional 
distribution of two particles separated by a distance n. 

Propostion 4. For two particles separated by a distance n = o(l/x), then 

• F[X(l,t + l),X(n-l,t + l)\X(0,t),X(n,t)] = \ + - + 0{x 2 ) 

4 7T 

• ¥[(X(-l,t + l),X(n + l,t + l)\X(0,t),X(n,t)] = \-- + 0{x 2 ) 

4 7T 

• F[(X(-1, t + l),X(n - l,t + l)\X(0,t), X(n,t)] = i + 0(x 2 ) 

• F[(X(l,t + l),X(n + l,t + l)\X(0,t),X(n,t)} = i + 0{x 2 ) 

Proof. The existence of particles at X(0,t) and X(n,t) means no dimers covering a edges at 
those locations. Thus, the dynamics of the particles are governed by the locations of dimers 
on the inverted triangle decoration. Therefore, the probability that there are particles at 
X(Q,t) and X(n,t) and they move closer together is given by, 

(5.16) 

( K^iv^va) K^_ n (v 4 ,v 4 (n,-n)) K~^_ n (v 4 , v 5 (n, -n)) \ 
K-\v 4 ,v 5 ) 

\ ••• o J 

where iT~^ n (v*, — n)) has been computed using LemmafTUl To compute the conditional 
probability that the particles move together given that there are two particles, it suffices to 
divide (I5.16|) by the joint probability of seeing two particles separated by distance n. This 
gives 1/4 + x/tt. The other quantities can be computed in a similar fashion. □ 

5.2. Location of particles in the scaling window. This subsection concentrates on 
finding the measure for the locations of particles along a horizontal line in the scaling window 
(x,x 2 ) for u = u c . In particular, we prove Theorem d] in the special case when u = u c or 
7 = "1- 

The proof of this theorem follows by a sequence of lemmas and propositions. Recall that 
Vi(n,m) is the vertex Vi in the (n,m) fundamental domain. In order to prove the above 
theorem, we need to find the entries of the inverse Kastelyn matrix in the re-scaling. This 
relies on the following: 

Lemma 11. Suppose that nx = a where a £ [0, oo). Then, 

(5.17) K'^vs, v 3 (n, -n)) = E x (a) + 0(x 2 ) 
Furthermore, 

(5.18) K- 1 (v3,v 4 (n,-n)) = E 2 (a)+E 3 (a)x + 0(x 2 ) 
where E\ (a) , E% (a) are defined in Theorem ^ and 

(5.19) E 3 (a) = — - 45/(2, 4a) + 45" L (-2, 4a) 

a7r 

where Bj(n, z) stands for the modified Bessel function of the first kind and Sl(ji, z) denotes 
the modified Struve Function. 
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The calculation is only given for K (vz,Vz(n, —n)) as the other term is analogous and 
can be completed with exactly the same steps. 

Proof. The first entry is given by 

1 f w^ 1 ^ n z n ^ 1 (-w 2 + z 2 ) 

(5.20) K" {v 3 ,v 3 (n,-n)) = 2 / ^-dzdw. 

(2-Kiy Jj2 (z - yi)(z - y 2 )(l + 2w- x l ) 

A residue can be taken with respect to z at y 2 due to |y2(^)| < 1 for all \w\ = 1. 
This results in an integrand containing the term + w) 2 + 4wx 2 in both the numer- 
ator and the denominator. To rewrite the integrand without this square root term, set 
x = 1/2^/— 2 — 1/r — r and take the change of variables w = (r — u 2 )/(l — u 2 r)). The 
transformation changes the contour of integration to a unit semi-circle from —1 to 1 pass- 
ing through the point— i. This contour shall be denoted C\. Note that r is the root of 
(1 + w) 2 + 4wx 2 whose absolute value is less than 1. Finally, setting R 2 = r gives 

r _ x / ^ 1 f 8R 2 (-1 + R 2 ) (R-u)(-l + Ru) 



(5.21) K- (t*,t*(n ,-»)) = — ^ - {R - - — } R _ 3R2 + R{3 + R2) u) f{w)dw 
where 

(R - u)(3R + R 3 - u - 3R 2 u) 



(5.22) f(w) ^_ 1 + Ru ^_ 1 _ 3R 2 + 3Ru + R 3 u ^J (_ U + R ( 3 + R 2 _ 3Ru yy 

The transformation, u = ZiXm > moves the contour C\ to a unit semi-circle from 1 to — 1 
passing through i. Denote this new contour by C 2 . Then, setting c x = — (4i? + 4i? 3 )/(l + 
6R 2 + i? 4 ) gives 

Tf _ l( , „ 8R 2 f {1 + R 2 - 2Rt) (t + R(-2 + Rt))(c x + tf' 1 

5.23 K (va, v^(n, — n)) = / j-i «t. 

V ^ V 3, 3V , 2?r . (l + 6 J R 2 + i? 4 ) 2 t"- 1 (c a; t + l) 1+n 

Setting n = a/x we can take a series expansion of the integrand in (|5.23p . This gives 
p m) ^ + xil + n . 2a . it + 2^ + 0(i2) 

Applying the bounded convergence theorem to (|5.23j) along with the change of variables 
t = e gives 

(5.25) «~- 1 (w 3 ,«3(n,-n)) = - / e " 4Qsine (1 + 2ixcos#(-l + 4asm6)) dO + 0(x 2 ) 

(5.26) =- e- 4asla6 d0 + O{x 2 ) 

TT 





The above integral can be evaluated in terms of the Bessel and Struve functions. 

□ 

The proof of Lemma [TT] also shows the 'correct' horizontal scaling. This can be shown 
by the series expansion of integrand in ()5.23j) . Taking a larger scaling, i.e. for e > set 
n = a/x 1+e , leads to (|5T2"3j) converging to zero. Taking a smaller scaling, i.e. set n = a/x e , 
then (|5.23p converges to 1 as x — > 0. In essence, n = a/x provides the right horizontal scaling 
to ensure decay of the inverse Kastelyn entry. 
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For 1 < i < m, let Xi S xTL C K denote possible particle locations on the lattice xZ. Let 
p™{x\, . . . ,x m ) denote the m point correlation function of seeing particles at (xi, . . . ,x m ) 
with the corresponding measure on R denoted by P" c , then 

Propostion 5. The m-point correlation function for the re-scaled particle locations, p m (yi, . . . . 
is given by 

(5.27) p* m (xi, ■ ■ ■ ,x m ) = Pf {M™( Xi , Xj ))™ =1 + 0(x 2m+2 ) 

where M^ es x is antisymmetric matrix with 



(5.28) 

and 
(5.29) 





-x + x 2 



-x - x 2 

7T 





xE 2 (\xj 



-xEi(\xj - xi\) 



-xE 2 (\xj - Xi\) - xE 3 (\xj - Xi\) 



X: 



+ xE : 



3{\Xj 



xEi(\xj 



Xi 



where Ei(a), E 2 (a) and E^(a) are defined in Lemma [771 

Proof of Proposition This is a direct consequence of Lemma[TT]and the underlying Pfaffian 
structure of computing the probability of the location of m particles on this particular dimer 
model. □ 

The above calculations provide a proof for Theorem [4] for the case u = u c . The proof 
follows an argument from |Bou07] . 

Proof of Theorem^ for u = u c . Notice that as and Pq c are both defined in the same 
space of measures, then it is enough to prove that P" c converges weakly to Pg c . Due to 
considering point processes, then tightness of the measures is guaranteed (see [DVJ88| ) . so it 
remains to show convergence of finite dimensional distributions. Let N{-) denote number of 
particles in a set, {A±, . . . , A^} denote a family of disjoint Borel sets in M and ni, . . . , G No, 
then we need to show 



(5.30) limP" c (iVMi) =m,...,N(A k ) 

x—tO 



n k ) = PS c (iV(Ai) = m, . . . , N(A k ) = n k ) 

in terms 



In order to prove such an equality, we need to write the measures P" c and n 
of their m-point correlation functions, which can be achieved by writing their respective 
generating functions. Denote n! = \\ 



Ek 

. . . z™* . The Taylor series of the generating function of IP 



nj, z = (zi, . . . , z k ) and z n 
" c about the point z = is 



(5.31) 



(5.32) 



Qx{z) = E a 



3=1 



"3 j 



pi 



n 



N(Aj)\ 



[(NiA-i-pj) 



Qx{z) can also be expressed as a power series with coefficients z n given by P" c (iV(Ai) 
rix, ■ ■ ■ , N(A k ) = n n ). The two representations of the generating function give 

-l) n d n 



(5.33) 



^(N(Ax) = nx,...,N(A k ) = n n ) 



nl 



dz" 



;Qx(z) 



«=(i,...,i) 
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Similarly, the generating function of P" c is given by 

k 



(5.34) 



Qo(z) = Yl 



P 



n 



N(AA] 



Therefore, (|5.30p is equivalent to showing 

d n {-\y 



(5.35) 



lim 



o dz n n! 



-Q x {z) 



L(iV(^)-Pi)! 
d n (-1) 1 



dz n n\ 



-Qo(z) 



z=(l,...,l) 



z=(l,...,l) 

Without loss of generality, suppose that the family of Borel sets, {Ai, . . . , A^} is open. 
Let A = {y £ 7, : yx £ A} for any set open set A. Using Proposition [5] gives 



(5.36) 



E, 



n 



iV(A,)! 



[(NiA^-nj)] 



J2 Pf(JW^(x 4 ,x i ))|" l ==1 



\n\ 



Taking the limit in x — > 0, the Riemann sums converge to 

(5-37) f ... f Pf (M™( m , yj ))W =1 dy? ...dyf 

which is exactly the m point correlation function for Pq c . By noting that all entries of 
(M^ x (xi, Xj))!™'^ are uniformly bounded by 1 (using Proposition[5]), we can use Hardamard's 
inequality, 



(5.38) 



E, 



n 



N(A,)l 



[(NiA^-nj)] 



<ni^i(2H) H/2 

3=1 



This implies that Q x (z) is an absolutely convergent series. For z in a compact set, then 
by the Lebesgue Dominated Convergence the derivatives of Q x (z) converge uniformly to 
Qo(z) in the limit x — > 0. This completes the proof of convergence of finite dimensional 
distributions. □ 



5.3. The Creation measure for u c . The creation singularities, under the scaling window 
(x, x 2 ) form a dense set of the plane. This can be seen by the fact that the expected number 
of creation singularities in a box of size (x -1 , x~ 2 ) is x" 1 . Here, we will only consider creation 
singularities which form non-trivial paths in the scaling window (x, x 2 ). By non-trivial paths, 
we mean paths of length e > in the scaling window (x,x 2 ). The measure on the set of 
creations, whose paths are non-trivial, is a locally finite (Lemma 1121) Poisson point process 
(Theorem [6]) . 

Let B be a box in (2xZ x 2x 2 Z) U ((2Z + l)x x (2Z + l)x 2 ) with corners (0, 0), (a, 0) and 
(0, (3) where < a < oo and < (3 < oo are fixed. Let H# denote the set consisting of all 
creation points in B for a realization of the underlying dimer covering. For (y, z) G let 
Pi(y,z,t — z) and P r (y,z,t — z) be the locations of the left and right paths emerging from 
(y, z) at time t > z. We use the convention that Pi(y, z) (or P r (y, z)) represents a realization 
of the left (or right) path starting from (y, z) £ Ss and \Pi(y, z)\ (or \P r (y,z)\) denotes the 
lifetime of the left (or right) path. Set, for a fixed e > 0, 



(5.39) 



B 



{(y,z) G E B : min(|P,(y, z)\, \P r (y, z)\) > e} , 
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This represents all the creations that create paths of length greater than e in the appropriate 
scaling window. The following shows that these creations are locally finite and do appear in 
the scaling window (x,x 2 ). 

Lemma 12. 

(5.40) < limElr^l < oo 

x— >0 

Proof. Notice that is less than or equal to the number of paths of length e. If lim^^o E 
oo, then this implies that the expected number of paths in the scaling limit is not locally 
finite, which is a contradiction. This proves lim^-s.o E | T B \ < oo. 

The inequality < lim-E^o E | T B | is more involved and requires the consideration of a 
random walk approximation. For (y,z) G T e B , let X t = P r (y,z,t — z) — Pi(y,z,t — z), with 
Xq = x and let Yt, with Y$ = 1, be the random walk which has iid increments of 1, 
and -1 with probability 1/4 — x/ir, 1/2 and 1/4 + x/ir. Let P r denote the corresponding 
probability measure for the random walk Yt started at r. Let = va.i{t > : Xt = 0} and 
T y = inf{t > : Y t = 0}. Then, 

Claim 1. There exists < T < Tq such that for < t < T , Xt stochastically dominates 
xY t/x 2. 

Note that T is the time of death for the shortest path emerging from [z, y) G T B and T > 
almost surely by definition of T e B . We are only interested in knowing that the difference of 
paths dominates a random walk for a strictly positive time. 

Proof of Claim. The evolution of each particle is a non Markovian process as each particle's 
trajectory is determined by the locations of the other particles. Proposition [5] shows that the 
pairwise force between two particles increases as the distance between the particles decreases. 
The candidate for the random walk approximation can be obtained from Proposition HI which 
is defined by Yt. Due to Yt having dynamics related to the maximum attraction between two 
particles, then the attraction between two paths P[(z,y) and P r (z,y) is always bounded by 
Yt, which does not matter on the locations of other paths. □ 

Define Y t (y, z) to be independent copies of the walk Y t , indexed by (y, z) £ Define 

(5.41) f B = {(y,z)eE B :T^ y ' z) >e/x 2 }. 

Due to paths from other creations can collide, then it is not immediately true that E[r^] > 
E[r^]. By Lemma \13\ the creation singularities on the 1/x x 1/x 2 scale are independent. 
Therefore, we only need to consider annihilations from paths with different creations that 
are not seen in the scaling window. For (yi, z±), (y2, £2) £ with y\ < yi and Z\ < z%, then 
Pi(y2, Z2, t — z\) — P r (yi, Zi, t — Zi) < —Yt stochastically. Therefore, we require 

(5.42) £P*/*(- Y */» = °) = P -i( y *A > ~ e l x \ Y - < °> Vr ) 

(5.43) = 0(x 2 ) 

which can be obtained by using Chebychev's inequality. Using this, we can bound the number 
of creations whose paths annihilate paths from other creations. This gives an O(l) bound, 
and hence, up to first order, E[r^] > E[r^] + O(x). This implies that 

(5.44) lim Efrfjl > lim E[f B ] 

x— >0 x— >0 
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and so the rest of the proof hinges on showing lim^^o E[Tb] > . 
By independence, 

(5.45) E |f B | > E[|H B |]Pi(T y > e/x 2 ) 
for some e > 0. Notice that 

(5.46) E[|3 B |] = l/x + 0(l). 

It remains to show that a lower bound for Pi(T y > e/x 2 ) is O(x). Using Kemperman's 
formula, then, 

(5.47) P!(T y >e/x 2 )= £ _L_ P (5 2r _ 1 = -1) 

r>e/x 2 

where Sp. = Yq H + Yjj. and 

(5.48) 



k=0 
2r- 1 



(5.49) 
Estimating gives 

[(2r-l)/4] +v ¥ / . w , ±lWl a 

,,50, r( S2r _ 1 = -i)> £ (^;i)( 2 \ +1 )Q) 2-(-> + o(^) 



(5.51) > - +0 — = \+0{x 2 ) 

~ V2^T \y/(2r-l)3) 

by performing a lower bound estimate on the smallest summand using Stirling's formula. 
Using the appropriate time re-scaling gives 

(5.52) Pi(T y > e/x 2 ) > C{e)x 

where C(e) represents some constant dependent on e. Substituting the required components 
back into ()5.45j) shows that lim^^o E[r B ] is nontrivial as required. 

□ 

Let B = lim x _ s> o-6 C M 2 . Let X denote the possible locations of creation singularities in 
B, so that for a dimer covering, X takes values in 3 B , Let Pi" c be the point process defined 
on B(M. 2 ), the Borel sets of M. 2 , such that the realization of all the creations whose paths 
are non-trivial in B is T e B . Here, the superscript V in the measure is to merely indicate 
singularity. The aim is to show weak convergence of the measure P|',e c to a Poisson point 
process. Before doing so, we first state and prove a lemma which was used in the proof of 
Lemma [12] and is used in the proof of Theorem (6) 

Lemma 13. For k > and yi, . . . ,y^ £ B, 



(5.53) P(X = yi , . . . ,X = y k ) = ]Jf(X = Vi ) + 0( 



x 2k+1 ) 



i=i 
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Proof. The proof heavily relies on the computation of the asymptotic inverse Kastelyn en- 
tries. These can be found algorithmically using a similar method employed in the proof of 
Lemma [TT1 The only difference is to take into account the vertical scaling. 

It remains to show independence up to a factor of x using the asymptotic inverse Kastelyn 
entries. Each creation consists of dimers covering edges v$ to i>i(l,0) and vq to ^(0, 1). 
Notice that 

(5.54) K- 1 (v 1 (0,l),v j ) = K- 1 {v 2 (l,0),v j ) + O{x) 

for j = 5, 6. For the off diagonal blocks, we have a similar result, namely 

(5.55) K~ l {v\ (0, l),Vj (n, m)) = K~\v 2 (l, 0), Vj (n, m)) + 0(x) 
for j = 5, 6, 

(5.56) # -1 («i(0, 1), v x (n + 1, m)) = K-\v*(\, 0), v x (n, m + 1)) + 0{x) 
and finally 

(5.57) # "^(O, 1), v 2 (n + 1, m)) = K- l (v 2 {l, 0), v 2 (n + 1, m)) + 0{x). 

This implies that the first two rows (and columns), apart from the 2 by 2 diagonal block, of 
the underlying inverse Kastelyn matrix ¥(X = y\, . . . , X = yj-), are the same up to 0(x). A 
similar property can be obtained for the third and fourth row. By considering a complete 
graph representation of the Pfaffian, see [God93| . gives 

(5.58) ¥(X = yi ,...,X = y k ) = (F(X = y x ) + 0(x 2 ))F(X = y 2 , . . . , X = y k ) 
Proceeding by induction gives the result. □ 



We can now state and prove a theorem about the locations of creations which give non- 
trivial paths in the scaling window (x,x 2 ). 



Theorem 6. For e > 0, Px,e c converges weakly to a Poisson point process in the limit x 
which has a finite non-zero intensity dependent on e. 



0. 



Proof of Theorem [H The proof follows the same recipe as the proof of Theorem 2] for u = 
u c (see Section 15. 2p . Again, tightness is already guaranteed and we only need to prove 
convergence of finite dimensional distributions. We borrow the multi-index notation from 
the proof of Theorem H] for u = u c . 

Let .Bi, . . . , B k represent disjoint boxes in (2xZ x 2x 2 Z) U ((2Z + l)x x (2Z + l)x 2 ) with 



Bi representing the box under the limit x — > for 1 < i < k. 
particle density for T e B . By Lemma \12\ we can write 



Let f m denote the m point 



(5.59) 



/ h(y)dy = lim V F(X = xi;n6 1^ 



x L eB 



limE 



u 



A 



< oo 



where P(-;xi G F € Bi ) is P(-) multiplied by the indicator xi G . Therefore, it is clear that 
F(X = x\\x\ G r^ i ) is 0(x 3 ). As the local dimer configurations are translation invariant 
(due to having zero boundary conditions), then J ^ fi(y)dy = /i(Bi)p where p is the intensity 
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of the point process. Using Lemma [131 and convergence of Riemann sums gives 
(5.60) 

N( B j)\ 



lim E^ c 

x^O 

(5.61) 

(5.62) 
Let 

(5.63) 

and set 
(5.64) 



n (N(B j 



n 3 )\ 



lim 



E 



"(X = xi;xi 6 T e Bl , . . . ,X = x\ n \;x\ n \ G T 



Bu 



xieBi,...,xi n t£B k 



lim 



E 



\x = x i; x! e r^)...p(jf = i w ;x|„| e rt )) 



xieSi^.^xi^ie-B,, 

h(y)dy 

Bi 



h(y)dy 



Qx(z) 



E 



Qo(*) = £ 

peNg 



peNg 
p! 



n 



B-, 



h(y)dy 



in 



h{y)dy 



i.e. a generating function for a Poisson point process. Notice that each generating function 
can also be written as a power series whose coefficients are given by the finite dimensional 
distribution probabilities for the corresponding measure on the sets B\ , . . . , B m . Following 
the same convergence argument given in the proof of Theorem |4] for u = u c , then it is clear 
that the derivatives of Q x {z) (with respect to z) converge to Qo(z). This proves the finite 
dimensional distributions convergence. 

□ 



6. The particle model for u 



In this section, we concentrate on the particle model for u = Ui and prove Theorem [3l The 
idea behind the proof is to form a bijection with the particle model and the color boundaries 
of a two color noisy voter model, defined in !( I M !)■">(. The weak convergence under the scaling 
window (x, x 2 ) is guaranteed by [FINR06J, who constructed the scaling window (x, x 2 ) of the 
noisy voter model. This limiting object is called the Continuum Noisy Voter Model. Finally, 
we show that the set of creations which give non-trivial paths in the scaling window (x, x 2 ) 
is a locally finite Poisson point process. Figure [9] shows a simulation of the particle model 
for x = 0.1 and u = Uj. 



6.1. Bijection with the Noisy Voter Model. In this subsection, we show that the par- 
ticles in the particle model have the same distribution as the boundaries of colors in two col- 
ored noisy voter model . This involves showing the distributions of particles are independent 
Bernoulli random variables, each trajectory is independent and creations are independent. 
Let 



(6.1) 



H(m, n) 



1 



(27Ti) i 



dz dw 



\w\=lJ\z\=l 



w n P(z,w) z w 
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Figure 10: A simulation of the particle model on a 100 by 100 grid with u = u\ and x = 0.1 
made using Glauber dynamics. 



Lemma 14. For m, n > 0, 

(-l) m+n (l - x) m (l + x) n {l + y/1 -x 2 ) 2+m ~ n {l -x 2 + Vl - x 2 )~ m+n 



3.2) H(m,n) 



Ax 3 



Proof. We only prove the result for m > n. To simplify calculations, set x to \/2r + 1/r. 
Computing the first integral via residue calculus gives 

1 f r(— 1 + 2r)w~ n (w — rw) r 



(6.3) / — : — ; = ttzoIw 

2-ki J\ w \ =1 2(-l + r + 2rw) m (-l + r + 2rw - w 2 + rw 2 ) 

The singularities of the integrand inside the region \w\ = 1 are located at w = (— r + 
\/2r — l)/(— 1 + r) and tt> = (1 — r)/(2r). The contribution from the residue at w = (— r + 
y/2r — l)/(— 1 + r) is given by 



+ v / 2^T\ m+n 



(6.4) i r v / 27^T ^ - j 

It remains to find the residue for (16. 3D at tu = (1 — r)/(2r). This can be achieved by 
splitting the term (—1 + r + 2ru; — w 2 + rw 2 ) = (r — l)(u> — ci)(u; — C2), where c\ = l/c% 
and |ci| < 1. Hence, the denominator in (16. 3[) can be separated into 

1 1/1 

.5) 



(— 1 + r + 2rw — w 2 + rw 2 ) 2\/2r — 1 \w — c\ W — C2. 
Each term can then be computed using the high dimensional residue formula: 

1 f rV2r -lw m - n (l-r) m 1 ,- ( r + ^/2r~ 



v ' 2iri J\ w+ ii=L\ =e 2(-l + r + 2rw) m {w - ct) 



2 2r 



and 



If rV2r - lw m - n (l - r) m 1 (r - 

' ' ' ~ 2~Ki J\ w+ ii^L\ =e 2(-l + r + 2rw) m (w - c 2 ) ~ 4 r 



2 2r 



m+n 
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Notice that contributions from (|6.4p , fj6.6[) and (|6.7|) sum up to give (|6.3p . As the contributions 
from (16.41) and (16.61) cancel out, then 



(6i 



1 



H(m,n) = — r\j2r — 1 



r - V2r - 1 
1 — r 



m+n 



Setting r = (1 + \/l — x 2 )/x 2 gives the result. □ 

Lemma[T4"lis an exact result. As a consequence, we can compute the entries of K~ 1 (vi, Vj(n, ■ 

exactly. In fact, it turns out that K~ l (v 3 , v 3 (n, — n)) = K~ l {y 3 ,V4 : {n, — n)) for u = 1 ~^J 2 ~ X ' 2 ■ 
This leads to an interpretation that particle locations are Bernoilli (-ffj), which is argued 
after the following proposition. 



Propostion 6. For n > 1, 

(6.9) 



F(X(0,t),X(n,t)) 



(1 + x) 2 

Proof. We can calculate the joint probability via the inclusion-exclusion formula given in 
Lemma [3l Hence, we require 

/ 1 - xK~ x (v 3 ,v^) -K- 1 (v 3 ,v 3 (n,-n)) -K' 1 ^, v 4 (n, -n)) \ 
— K _1 (u4, vs(n, — n)) — i^ _1 (f4, v±{n, — n)) 

l-xK- x (v z ,v A ) 





(6.10) 



Pf 



V 



/ 



where Vi for i G {1, ... ,6} represents a vertex i in the fundamental domain and Vi(m,n) 
represents vertex i a distance (m, n) from the fundamental domain. 

Using residue calculations, K~ l {v^,Vi) = l/(x(x + 1)). It remains to calculate the in- 
teraction between the two edges, which can be found by invoking Lemma [T4l Due to 

, by symmetry (16,101) becomes 

l-xK-\vs,v$ 0) \ 





K l (v 3 ,v 3 (n, -n)) = -K 1 (v 3 , u 4 (n, -n)) for u = 
/ 



(6.11) 



Pf 



V 



0) 

K^Vi, v 3 (n, -n)) K' 1 ^, v 3 (n, -n)) 
1 -xK^iy^vt) 





The Pfaffian of the above matrix represents counting the number of weighted perfect match- 
ings of a graph with four vertices w\, . . . ,W4, with the edges weighted by the above matrix. 
Clearly, we can only have one possible contribution (the edges W1W2 and u^w^), which gives 
the desired result. □ 

The above method can be generalized to consider the locations of an arbitrary number of 
particles. 

Propostion 7. For %\ < • • • < i n , 

x n 

(6.12) ¥(X(h, t), . . . , X(i n ,t)) = Tf^p 

Proof. Using the inclusion-exclusion formula given in LemmaEl we can write W(X(i\, t), . . . , X{% 
as a Pfaffian of a matrix with block diagonals given by 

1 - xK- l (v 3 ,Vi) 



3.13) 
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.14) 



and off-diagonal blocks given by 

-K~ l (v 3 ,v 3 {i k , -ik)) -K- 1 (v 3 ,v i (i k , -ik)) 
-K- l {v A ,v 3 (i k , -ik)) -K~ l {v^Vi(ik, -ik)) 

for k € {1, . . . , n}. Using the fact that —K~ 1 (v 3n v 3 {ik, —ik)) = ^ (i>4, vz{ihi —ik)), then by 
applying the graphical representation of the Pfaffian inductively (described in Proposition 
[6]), we can write 

(6.15) nX(h,t), . . . , X(i n ,t)) = (1 - xK~\v 3 , Vi )) n . 

which gives the required result. □ 

The locations of the particles are independent and Bernoulli (■jtj). It now remains to look 
at the dynamics. 



P(L ,L n ) =«VPf 



Propostion 8. For an arbitrary horizontal line, dimers over the b edges are iid with distri- 
bution Bernoulli ({1/2 — ^z?L ). 

Proof. We consider an arbitrary horizontal line. Suppose that Li represents the event that 
a dimer covers a b edge going left (i.e. from ^5 to v\) where i represents the position of V4 
on the horizontal axis. Let Ri denote the event that a dimer covers a b edge going right (i.e. 
from vq to vi) where i represents the position on the lattice V4 on the horizontal axis. 

We first consider the probability of two left going b edges a distance n away. Without 
loss of generality, we can consider one of these edges to be located at the origin. 

(6.16) 

/ K- 1 (v 5 ,v 2 (l,0)) K- l {v 5 ,v 5 (n,-n)) R- l (v 5 , v 2 {l - n, n)) \ 

K~ l (y 2 ,v 5 ) K-\v 2 ,v 5 ) 

K- 1 (v 5 ,v 2 {0,l)) 

V ••• J 
where we have used '. . . ' to denote the lower triangle entries (the matrix is anti-symmetric). 
Then, explicit calculation using Lemma [14] gives 

(6.17) K-\v 5i v 5 (n, -n)) = K~ l {v^,v 2 {n, 1 - n)) 
Therefore, using row and column operations on the above matrix gives 

F(L ,L n ) = u 2 x 2 (K }(5,2)) 2 

1-Vl-x 2 
2(1 + x) 

A similar result for F(Ro,R n ) is true, but is reliant on 

(6.18) K-\v & ,v & {n,-n)) = K~\v &1 v 1 {n + l,-n)). 

The above calculation can be generalized to find the joint probabilities of observing dimers 
covering any collection of rightward and leftward b edges located on the same horizontal 
line. Again, this joint probability is just the product of the probabilities of observing a 
dimer covering each edge. Furthermore, we can show that the probability of observing n 
creations along a horizontal line, is given by the product of the probabilities of dimer covering 
each edge. For example, 

(6.19) P(L , Ro, L n , Rn) = i£x\K- x (vs, v 2 (0, l))) 2 ^ 1 ^, ^(0, l))) 2 
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□ 



A further calculation shows that the trajectories of particles are independent. This relies 
on computing the possible trajectories of the particles given their locations. This calculation 
follows in the same vein as the previous propositions. 

Proof of Theorem For u = 1 ~^~ rr2 , the dimer covering of the Fisher lattice can be 
thought as follows: 

• Take an arbitrary horizontal line along the a edges. With product measure, we can 
place dimers on a edges with an iid Bernouilli (ytpj) distribution. 

• If there is no dimer covering an a edge on the row below (or above), then choose a 
dimer to cover either the left or right b edge directly above (or below) with probability 
^ . This is a direct consequence of seeing one b edge conditioned on there not being a 
b edge on the row below (above). The surrounding dimers are automatically selected 
from this probabilistic choice. 

• If there is a dimer covering a edge on the row below, then with probability (1 — 
y/1 — x 2 )/(2 + 2x), select two dimers to cover both directly above (below) b edges. 
Otherwise, cover the edge v$ to Vq with a dimer. 

As a consequence of locating all the dimers on the b edge row, determines the existence of 
dimer on the above (below ) a edges row. 

Due to the independence of the entries, the model can be viewed as a one dimensional 
noisy voter model with two colors, say red and blue. The particles in the dimer model 
represent the boundary between the two colors. A site contained within a connected region 
of the same color can flip color with Bernoulli(l — \/T— x 2 ) / (2 + 2x). This represents the 
noise. Hence, the particle model is in direct correspondence with the noisy voter model. The 
scaling window (x, x 2 ) of the noisy voter model was constructed in [FINR06] and named the 
Continuum Noisy Voter model. 

□ 

An equivalent interpretation of this dimer model is a one-dimensional Ising Model with 
heat bath dynamics. This can be seen via the two-dimensional Ising model representation 
of the dimer model. The fact that there are no correlations amongst the spins along any 
horizontal line of the Ising model in two dimensions implies that the interaction of the spins 
is only dependent on the previous levels. 



6.2. Creation Measure. Similar to u = u c , the set of creations for the case u = Uj forms 
a dense set in the scaling window (x,x 2 ). However, we can restrict to the set of creations 
which give paths of length greater than e in the scaling window. This section shows that 
this random set is finite and its points form a Poisson point process. 
We use the same notation as Section [5.31 without further mention. 

Lemma 15. 

(6.20) < limElr^l < oo 

x— >0 

Proof. The proof of the upper bound follows by the same argument given in Lemma [T2l 

For the lower bound, define X t (y, z) = P r (y, z,t — z) — Pi(y, z,t — z) for (y, z) G Sg. This 
means {Xt(y, z) : (y, z) 6 Sg} is a set of independent walks with i.i.d. increments of 1 with 



40 SUNIL CHHITA 



probability 1/4, with probability 1/2 and -1 with probability 1/4. Let Xt represent one 
copy of these walks. Then, using the same arguments of Lemma \12\ 

(6.21) lim E[|r^|] = lim E[|S B |]Pi(r^ > e/x 2 ) 

x— >0 x— s>0 

¥\{Tq > e/x 2 ) has already been computed in Lemma PT2l and is Q(x). E[|Hb|] is 1/x 3 mul- 
tiplied by the expected number of creations in a fundamental domain. Therefore, E[|E#|] = 
c/x + 0(1). Hence, 

(6.22) Urn i E[\T%\] > 

x— H) 

□ 

As in Section E31 E[|r^|] is dependent on e. Moreover, lim e _>o E[|T^|] = oo. 

Theorem 7. P^' e * converges weakly to a Poisson point process as x — > 0. 

Proof. By using Lemma [T4"l then it can be shown that the creation singularities are inde- 
pendent. In other words, the statement Lemma [13] can be replaced by 

k 

(6.23) P(X = y u ...,X = y k ) = l\¥(X = Vi ) 

i=l 

for k > and yi, . . . ,yk G B. Following the proof of Theorem [6j using the above equation 
where necessary gives the required result. □ 

Remark: Let Pf = {Pi(z,y) : (z,y) G T e B } and define P r e accordingly. Let P* 1 ™ 9 de- 
note the limiting Poisson point process for Px,e 9 with e > fixed. Then, define R(e) = 
{Fq £ 9 ,Pi,P^} denotes the CNVM conditioned on having path lengths of at least e > 0. 

7. The scaling around u r and m 



This section concentrates on the behavior exhibited when u = (1— \J\ — 7X 2 ) / '(7a: 2 ) with 
7 G K\{— 1,1} and I7I = o(l/x)(i.e. u = 1/2 — jx 2 /8). We find the stationary measure 
for the locations of particles on a horizontal line for the scaling window (x,x 2 ). There are 
three cases to be considered, namely, 7 < —1, 7 G (—1, 1) and 7 > 1. These correspond to 
u < u c , u c < u < Ui and u > U{. For each case, particle locations are given by a Pfaffian 
point process. However, each process exhibits different local behaviors, which are discussed 
after the following lemma. 

Lemma 16. Suppose that nx = a where a G [0, 00) and 7 G — 1, 1}. Then 

(7.1) K-\v3, v 3 (n, -n)) = Ej(a) + 0(x) 
and 

(7.2) K~ l {vz, v A (n, -n)) = E%(a) + 0(x) 
where 

1 f e i(7) 9 -ap 

(7.3) Ej(a) = -— P ^dp + Q{x) 

We 2 ( 7 ) ^4 + 87 + 4 7 2 - 12p 2 + 4 7 p 2 + 

and 

1 r e ^(~l) (—0 — 2-v — n 2 )p~ a P 

( 7 - 4 ) ^» = ~ / n } = = = 4 I o J P + 0(x) 



TTl 



<-2 



( 7 ) 2^4 + 87 + 4 7 2 - I2p 2 + 4 7 p 2 + p 4 
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where ei( 7 ) = 2 + ^2(1 - 7 ) and e 2 ( 7 ) = (2 - ^2(1 - 7 ))I 7 >-i + (-2 + ^2(1 - 7 ))I 7 <_i . 

Lemma [16] provides the inverse Kastelyn entries, which can be used to compute cor- 
relations. For 7 > 1, notice that ei( 7 ) = e2( 7 ). For 7 < 1, it can be shown that 
|i?^(a)| > I-Evjf^)!- This implies that any pair of particles is positively correlated. How- 
ever, for 7 > 1, it can be shown that |i?^(a)| < l-E^aOli which shows that any pair of 
particles is negatively correlated. 

Lemma [16] is also valid for 7 = — 1 but not for 7 = 1. From the observation that 
62 ( — 1) = 0, then clearly the decay of E^ 1 (a) in the limit a — > is slower than the decay 
of Ej(a) in the limit a — > for other 7 { — 1, 1}- This results in the covariance to decay 
slower in distance at u c , then any other value of u. 

The entry for K _1 (v3, v^) can also be computed up to O(x) in a similar fashion as the 
other entries in Lemma [TBI . Let e( 7 ) = lim^^o — X K~ 1 ( V 3, V 4 : ))- Then, 

1 , 1 f e i(T) 9 — 2-v -I- n 2 

(7.5) e( 7 ) = lim - (1 - xK~\v^v,)) = - / - 7 +P dp 



£ ' 7ri 7 e2 ( 7) 2^4 + 8 7 + 4 7 2 - 12p 2 + 4 7 p 2 + p 4 

Under the scaling window (x,x 2 ) and the limit x — > 0, the particle locations are given by 
Theorem HI We can now prove Theorem O 

Proof of Theorem^ The structure of the matrix comes from Lemma [3] The relevant entries 
come from Lemma [TBI The proof of weak convergence is the same as the proof in Theorem 
U]for u = u c . □ 

It remains to compute the inverse Kastelyn entries. Let 

(7.6) P(z,w) = zwP(z,w) = zwdetK(z,w) 

Proof of Lemma\Tb\ For the sake of brevity, we only compute K -1 ^, v 3 (n, —n)). The entry 
is given by 

1 y _ u 2( u; -n+l n-l _ ^-n-l^n+l^ 



(7.7) K L (v 3 ,v 3 (n,-n)) = -—2 * g- ^dzdn; 

(27TZ) 2 J T 2 P{Z,W) 

Taking the transformation z = C w gives 

(7.8) K -i (v3Mn ^ n))= VLk ^ ^ rfwdC , 

(27T«)^ J T 2 P((w,w) 

As P((w,w) = C(w — wi)(w — W2) where C, w;i and tU2 are rational functions of £, with 
|wi| > 1 and \w2\ < 1 then computing the integral with respect to w via residue calculus 
gives 

(7-9) K-\v 3 Mn-n)) = [ ^ ^^L' 1]) d( 



2m y K | =1 ^^(C) 

where 

(7.10) A(() = (wi - w 2 )u((l + C)(-l + u 2 x 2 ) 



Recalling the expressions ei( 7 ) and 62(7) and denote 63(7) = 3 — 7 + 2-^/2(1 — 7) and 
64(7) = 3 — 7 — 2^/2(1 — 7). Then, A(Q can be factorized into the following 
(7.11) 

A{() = u 4 (-l+x 2 ) 2 (C-(l+eix+e3X 2 ))(C-(l+e2X+e 4 a; 2 ))(C-(l-eix+e3X 2 ))(C-(l-e2X+e4X 2 ))+0(x 3 
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Taking the change of variables £ = 1 — x£ transforms -A(C) in (|7.1ip into the following 
(7.12) 

A(0 = x 4 u 4 (-l+x 2 ) 2 ((£ - (ei + e 3 x))(£ - (e 2 + e 4 s))(£ - ( ei + e 3 x))(£ - (e 2 + e A x)) + 0(x 2 )) 
while (17.91) becomes 



(7.13) iT 1 ^ 1*3(71, -n)) = — / ^ ^ " ~ 



2vri 7^-1/^1=1/^ \M(0 

The contour of integration can be pushed to a contour independent of x surrounding the two 
roots of -A(£) with positive real part, namely a contour surrounding ei(-y) and e 2 (7). Notice 
that ei(7) and e 2 (7) are defined so that we can keep track of the roots with positive real 
part. Then, the numerator of the above integrand tends to 

(7.14) ^e~ a ^x 2 
as x —?• whereas the denominator tends to 

(7.15) - -x 2 ^/4 + 87 + 4 7 2 - 12£ 2 + 4 7 £ 2 + £ 4 
By the denominated convergence theorem, we obtain 

(7.16) K- 1 (v 3 ,v 3 (n,-n)) = / — z d£. 

V ' 2m J ^4 + 87 + 4 7 2 - 12£ 2 + 4 7 £ 2 + £ 4 

where the integral is over a contour surrounding 61(7) and e 2 (7). Splitting up the contour 
into two line integrals gives the required result. 

□ 

Although Lemma [TBI incorporates Lemma [TT1 the technique used to calculate Lemma [TT] 
is more general. Moreover, we can use the same technique found in Lemma [TT] when the 
indices of z and w are different. This is required in Lemma [T3l The calculation for Lemma 
[TBI cannot be generalized in this way. 

7.1. Correlation Length. This subsection focuses on proving Lemma 11.161 Recall that 
the correlation length, £(7) ( II . 16f) is defined in terms C(a,j) = Ej(a) 2 — E^ia) 2 for 7 7^ 1 
and C(a, 1) = for all a > 0. 

Proof of Lemma \1.16l We split the proof into the various cases. 

For 7 = 1, C(a, 1) = for all a > 0, which implies that £(1) = 0. As C(a, -1) ~ 1/a 2 , 
then £(— 1) = 00. 

For — 1 < 7 < 1, we can set up bounds for El (pi) and E^ipc) for a large: 

(7.17) e- aei ^Ej(0) < Ej(a) < e - ae ^ Ej(0) 
and similarly 

(7.18) e- aei W^(0) < E%(a) < e^^-E^O) 
Therefore, 

(7.19) e~ 2aei ^Ej(0) 2 - e- 2ae2 ^E2(0) 2 < C(a, 7 ) < e~ 2ae ^ Ej{Q) 2 - e~ 2aei ( ^E]{0) 2 . 
Using these bounds in the definition of £(7) gives 

(7.20) — \-r < £(7) < 



2e 2 (7) 2e 2 (7) 
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as required. 

For 7 < — 1, (|7.17p still holds but (|7.18p does not necessarily hold. Split E^a) into its 
positive and negative parts, i.e. 

(7.21) £ 2 >) = (£ 2 >)) + -(£ 2 >))_, 



with (.Eg (a))+ having bounds of integration 62(7) up to y/2(--y - 1) and (£^(a))_ having 
bounds of integration yj2(— 7 — 1) up to 61(7). Then, 

(7.22) \E<{a)\ < e-^ a {E<m + . 
and 

(7.23) \El{a)\ > e-V 2 (-T-i)°(^(o)) + + e" ei W°(^(a))_ 

Plugging these bounds into C(a, 7) and using the definition of <!;(«, 7) gives the required 
result. 

For 7 > 1, the underlying integrals are contour integrals, with the contour having constant 
real part. From the definition of C(a,j), we can consider \Ej(a) — Ej(a)\ and \Ej(a) — 
i£j(a)| separately. For notational convenience, set r = \/2( / y — 1). By applying the change 
of variables p 1— > 2 + then 



(7.24) I #7 (a) - El{a)\ = 2— Re /"" W- (p 0(p + r) 

V ' ' U 7 2K n 7T 7o y (p-2i + r)(p-2i + r) 

and 



(7.25) \Ej(a) - Ej(a)\ = 2^Re f / (P ~ ^ + - 2^ + r) ^ 

1 2 y_ r y (p_ r )(p +r ) 

by splitting up the domain of integration into two parts (— r, 0) and (0, r). The integral in 
(|7.24p decays at rate 1/r. This can be seen by applying integration by parts and integrating 
the term y/p — re~ tap . Employing a similar technique shows that (|7.25p decays at rate 1/ ' \fr. 
Therefore, 

(7.26) lim - log C(a, 7) = -\. 

Plugging back into the definition of the correlation length gives the result. 

□ 



We conclude with some remarks about the correlation length. Choosing 7 = 1 + Cx a , 
where C > and a > 0, forces the covariance between particles to converge to in the scaling 
window (x,x 2 ). This can be shown either using an analogous calculation (|7.17p . or by using 
an expansion of the covariance around u = U{. As lim T _ >1 + £(7) = 1/4, lim 7 _ s>1 - £(7) = 1/4 
and £(1) = 0, then there is a point of discontinuity for the correlation length at 7 = 1. 
This is a consequence of the definition of the correlation length as opposed to a difference 
of measures. This can be seen by setting u = Ui + ex 2 and showing that all the finite 
distributions converge (after taking the scaling window) in the limit e — > 0. 
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8. Concluding Remarks 

Here, we describe some possibilities for future work that we have been unable to complete 
so far. 

As noted before, we have been unable to find the limiting measure for the scaling window 
(x, x 2 ) when u = (1— yl - 7X 2 ) / (jx 2 ) for 7 7^ 1. When 7 = — 1, we expect that this limiting 
measure is scale invariant but not rotationally invariant. This is due to the underlying dimer 
model (or Ising model) being anisotropic. The techniques used in this paper have primarily 
concentrated on spin-spin correlations. In the same way that the scaling limit of percolation 
(see |CN06| ) required SLE and CLE (see [Wer04| or [LawOSj for fine surveys of SLE), we 
expect that the scaling windows of this Ising Model requires some family of random curves 
and loops. 

Taking x = 1 in the underlying dimer model forces the coefficient of z/w and w/z in 
the characteristic polynomial to be zero. For u = 1, P(z,w) is a constant. This implies 
that spins from the Ising Model are chosen independently with probability 1/2, which is 
exactly percolation on the triangular lattice. For u < 1, the inverse Kastelyn entries can be 
written in terms of the Green's function for random walks on the triangular lattice. Choosing 
u = \/2 — 1, the Green's function is the inverse of the Laplacian. 

If we take u = Ui + e 2 where e is order x a for a S (0, 1), then 

(8.1) E[N b ] = E[N a ] + 0(x) = 0(e) 
and 

(8.2) E[N X ] = x 2 + eO(x 2 \ogx). 

These can be shown by expanding out the appropriate elliptic integrals defined in Section 
[3l By taking the scaling window (x a , x 2a ) implies that the density of particles is of constant 
order. For a < 2/3, there is a dense set of creations in the scaling window. For a > 2/3, 
then there are no creations or annihilations seen in the scaling window (x a , x 2a ). We expect 
the limiting measure to be Dyson Brownian motions. For a = 2/3, the creations have a 
finite density. Roughly speaking, we expect the repulsion of the particles to be weak enough 
to allow some annihilations. 

Finally, is there a random matrix ensemble, whose eigenvalues have the same distribution 
as the particle locations in the scaling limit for u = u c ? This question is motivated from the 
fact that dimer coverings on the honeycomb lattice are a 'discretized' version of a GUE. 
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